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1 .  Stress  Measurement  In  Single  and  Polycrystalline  Ceramics  using  Their  Optical 
Fluorescence,  Q.  Ma  and  D.  R.  Clarke,  Journal  of  the  American  Ceramic  Society, 
submitted. 

2.  Direct  Measurement  of  Stresses  in  Sapphire  Fibers  In  Composites,  Q.  Ma  and  D.  R. 
Clarke,  Acta  Metaliurgica  et  Materialia,  to  be  submitted. 

3.  Optical  Fluorescence  From  Chromium  Solutes  In  Sapphire:  A  Probe  of  the  Image 
Stress.  Q.  Ma  and  D.  R.  Clarke,  Acta  Metaliurgica  et  Materialia,  to  be  submitted. 

4.  Acoustic  Emission  During  Polarization  Switching  in  Ferroelectric  Ceramics,  W.  Pan 
and  H.  C.  Cao,  Ferroelectrics,  Vol.  129,  pp.  119-126  (1992). 

5.  The  Structure  and  Properties  of  a  Ferroelectric  Electrode  Interface,  H.  C.  Cao,  M. 

De  Graef  and  A.  G.  Evans,  submitted  to  Journal  of  the  American  Ceramic  Society. 

6.  Electric  Field-Induced  Fatigue  Crack  Extension  in  Ferroelectric  Ceramics,  H.  C.  Cao 
and  A.  G.  Evans,  to  be  completed. 

7.  Singular  Behavior  Near  a  Crack  Tip  in  an  Elastic  Electrostrictive  Material,  by 
H.  C.  Cao  and  R.  M.  McMeeking,  to  be  published. 

8.  Non-Linear  Constitutive  Properties  of  Piezoelectric  Ceramics,  H.  C.  Cao  and  A.  G. 
Evans,  submitted  to  the  Proceedings  of  the  MRS  Spring  Meeting . 

9.  Non-Linear  Deformation  of  Ferroelectric  Ceramics,  H.  C.  Cao  and  A.  G.  Evans, 
submitted  to  Journal  of  the  American  Ceramic  Society 
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DIRECT  MEASUREMENT  OF  STRESSES  IN  SAPPHIRE  FIBERS  IN 

COMPOSITES 

Qing  Ma  and  David  R.  Clarke 
Materials  Department 
University  of  California 
Santa  Barbara,  California  93106 

Abstract 

The  residual  stresses  in  sapphire  fibers  embedded  in  a  number  of  ceramic  and  inter- 
metallic  matrices  are  measured  non-destructively  and  in-situ  using  the  stress- 
induced  frequency  shift  in  fluorescence  from  chromium  dopants  in  the  fibers.  The 
stresses  as  a  function  of  depth  into  the  fibers  from  the  free  surface  are  measured  and 
from  the  depth  profile  we  have  determined  the  interfacial  shear  stresses  between 
the  fibers  and  the  matrices.  The  stress  distributions  provide  a  validation  of  the  shear 
lag  mechanics  models  of  Hutchinson  and  Jensen,  as  well  as  a  test  of  finite  element 
computations. 


SIKESS  MEASUREMENT  IN  SINGLE  A^JD  POLYCRYSTALLINE  CERAMICS 
USING  THEIR  OPTICAL  FLUORESCENCE 


Qing  Ma  and  David  R.  Clarke 
Materials  Department 
University  of  California 
Santa  Barbara,  California  93106 


A  general  methodology  is  developed  for  determining  the  state  of  stress  and  the 
numerical  value  of  the  stresses  from  observed  shifts  and  broadening  of  optical 
fluorescence  lines.  The  method  is  based  on  the  piezo-spectroscopic  properties  of 
single  crystals.  We  present  general  relationships  between  the  measured  fluorescence 
shifts  and  the  stress  state  for  a  number  of  illustrative  cases,  pertinent  to  both 
polycrystalline  and  single  crystal  ceramics  under  stress.  These  include  measuring  the 
stresses  applied  to  poycrystalline  ceramics,  the  residual  stress  distribution  due  to 
crystallographically  anisotropic  thermal  expansion,  and  the  stresses  applied  to  single 
crystals.  Using  the  recently  implemented  technique  of  performing  the  fluorescence 
measurements  in  an  optical  microprobe,  we  also  provide  experimental  tests  of  the 
relationships  derived. 


Submitted  to  the  Journal  of  the  American  Ceramic  Society 


optical  Fluorescence  From  Chromium  Ions  In  Sapphire: 

A  Probe  Of  The  Image  Stress. 

Qing  Ma  and  David  R.  Clarke 

Materials  Department, 

University  of  California, 

Santa  Barbara,  CA  93106 

Abstract 

The  frequency  of  fluorescence  from  chromium  substitutional  solute  in  sapphire  is 
shown  to  be  a  direct  probe  of  the  image  stress.  This  provides  a  confirmation  of  the 
existence  of  the  image  stress,  a  component  of  the  total  internal  stress  due  to 
inclusions,  twins  and  precipitates,  first  introduced  by  Eshelby.  The  image  stress 
contribution  can  be  exploited  in  using  chromium  fluorescence  as  a  sensor  of  stresses 
in  solids. 


To  be  submitted  to  Acta  Metallurgica  et  Materialia 
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STRESS  \IEASUREMENT  IN  SINGLE  AND  POLYCRYSTALLINE  CERAMICS 
USING  THEIR  OPTICAL  FLUORESCENCE 


Qing  Ma  and  David  R.  Clarke 

Materials  Department 
University  of  California 
Santa  Barbara,  California  93106 

Abstract 

A  general  methodology  is  developed  for  determining  the  state  of  stress  and  the 
numerical  value  of  the  stresses  from  observed  shifts  and  broadening  of  optical 
fluorescence  lines.  The  method  is  based  on  the  piezo-spectroscopic  properties  of 
single  crystals.  We  present  general  relationships  between  the  measured  fluorescence 
shifts  and  the  stress  state  for  a  number  of  illustrative  cases,  pertinent  to  both 
polycrystalline  and  single  crystal  ceramics  under  stress.  These  include  measuring  the 
stresses  applied  to  poycrystalline  ceramics,  the  residual  stress  distribution  due  to 
crystallographically  anisotropic  thermal  expansion,  and  the  stresses  applied  to  single 
crystals.  Using  the  recently  implemented  technique  of  performing  the  fluorescence 
measurements  in  an  optical  microprobe,  we  also  provide  experimental  tests  of  the 
relationships  derived. 


1  INTRODUCTION 
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Despite  the  importance  of  measuring  the  stress  state  in  a  materiad  or  structure  there  are 
few  quantitative  methods  available.  Those  that  have  been  developed  and  widely  used,  such  as 
photoelastic  methods  and  X-ray  methods,  whilst  accurate  are  generally  limited  to  examining 
flat  surfaces  and  ones  where  the  area  interrogated  is  not  obscured  or  shadowed  by  features 
such  as  steps.  In  addition  they  axe  difficult  to  apply  to  the  measurement  of  stresses  from 
small  areas,  a  few  microns  across,  and  are  not  generadly  microscope-based  thus  precluding 
making  measurements  from  featu'es  identifiable  during  examination  of  a  microstructure. 

For  these,  and  related,  reasons  a  number  of  opticcil  techniques  based  on  piezo-spectroscopic 
phenomena  have  been  proposed.  One  gaining  increased  support,  especially  in  VLSI  silicon 
based  technologies,  is  that  of  Raman  spectroscopy  using  an  optical  microprobe  [1].  The 
Raman  lines  of  silicon  at  520  cm~^  axe  shifted  by  the  application  of  a  stress,  and  so  re¬ 
gions  on  a  silicon  wafer  around  fine  metal  lines  and  dielectric  features  can  be  probed,  with 
approximately  a  micron  resolution,  and  their  stress  evaluated.  Another  piezo-spectroscopic 
technique  proposed  in  the  late  1970’s  for  the  me2isurement  of  stress,  is  that  of  fluorescence 
spectroscopy  [2].  The  suggestion,  made  by  Grabner,  followed  the  demonstration  at  the  Bu¬ 
reau  of  Standards  that  the  pressure  in  high-pressure  experimental  cells  could  be  calibrated 
by  the  shift  of  the  chromium  fluorescence  lines  from  fine  particles  of  ruby  embedded  in  the 
high-pressure  cell  [3].  In  turn,  this  work  built  on  the  luminescence  studies  of  chromium 
doped  sapphire  by  Schawlow  [4]  and,  independently  by  Sugano  and  Tanabe  [5],  as  part  of 
their  investigation  of  possible  Iciser  materieds.  They  found  that  the  luminescence  lines  in 
chromium-doped  sapphire  were  shifted  by  an  applied  stress.  We  have  recently  demonstrated 
that  by  implementing  the  fluorescence  spectroscopy  using  an  optical  microprobe  it  is  pos- 
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sible  to  measure  stresses  in  chromium-doped  alumina  from  regions  identified  in  the  optical 
microscope  and  with  a  spatial  resolution  of  a  few  microns  [6]. 

In  this  paper  we  present  a  method  of  analyzing  an  observed  frequency  shift  in  terms  of 
the  piezo-spectroscopic  coefficients  and  the  stress  (and  hence  strain)  on  the  crystal  structure 
containing  a  fluorescing  atom.  This  provides  a  correction  to  the  equation  derived  by  Grab- 
ner  for  the  case  of  chromium  fluorescence  in  2dumina  [2].  We  then  extend  the  einaiysis  to 
the  measurement  of  stresses  in  inhomogeneously  stressed  materials,  presenting  a  number  of 
illustrative  examples  of  the  application  of  the  fluorescence  technique  to  the  measurement  of 
commonly  occurring  types  of  stresses  in  ceramics.  These  provide  the  basis  for  future  work  in 
the  widespread  application  of  the  technique  to  measure  stresses  in  ceramic  composites  and 
at  ceramic -interfaces. 

2  PRINCIPLES  OF  PIEZO-SPECTROSCOPIC  ANALYSIS 

The  basis  of  all  piezo-spectroscopic  methods  for  the  measurement  of  stress  in  crystalline 
materials  is  that  an  applied  stress  strains  the  lattice  and  alters  the  energy  of  transitions 
between  electronic  or  vibrational  states.  In  this  way,  the  energy  of  any  radiative  transitions, 
such  as  those  that  give  rise  to  luminescence,  also  vaxies  and  produces  systematic  shifts 
in  the  frequency  of  luminescence  lines  such  as  fluorescence  lines.  In  common  with  other 
spectroscopies,  the  sharpness  of  the  fluorescence  line  depends  on  temperature  emd  variations 
in  the  site  position  of  the  fluorescing  atomic  species.  When  the  fluorescing  ion  is  located  on 
a  lattice  site  in  a  tightly  bound,  high  elastic  modulus  material,  £is  is  the  case  for  chromium 
substitution  for  an  aluminum  ion  in  the  sapphire  lattice,  the  lines  can  be  exceptionally  sharp, 
broadened  only  by  phonon  scattering.  Interstitial  dopants,  where  there  is  greater  variability 
in  the  dopant  site  or  in  a  random  network  structure,  such  as  chromium  in  glasses,  give  rise 
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to  much  broader,  and  usually  weaJcer,  luminescence  lines. 

The  principle  of  relating  an  observed  line  shift  in  a  fluorescence  or  absorption  spectrum 
to  the  state  of  stress  have  been  described  previously  by  Grabner  [2]  in  reference  to  the 
measurement  of  stresses  in  sintered,  polycrystalline  alumina  containing  a  trace  of  Cr^'*’.  The 
more  general  case  is  presented  here. 

When  a  crystal  is  subject  to  a  homogeneous  stress  (in  crystallographic  frame  of  refer¬ 
ence),  the  change  in  frequency,  Ai/,  of  a  luminescence  line  is  given  by  the  tenscrial  relation¬ 
ship: 

A:/  =  n.j(T*  (1) 

where  Ilij’s  are  the  piezo- spectroscopic  coefficients  relating  frequency  to  stress.  They  form 
symmetrical  second  rank  tensors,  which  have  symmetry  properties  governed  by  the  point 
symmetry  group  of  the  fluorescing  ion.  They  may  also  have  different  numerical  values  for 
different  fluorescence  line  energies.  No  appreciable  stress  induced  broadening  heis  been  re¬ 
ported  at  room  temperature  so  the  lines  merely  shift  with  stress  without  changing  shape. 
(Stress  induced  splitting  is  detectable  at  liquid  helium  temperatures  but  is  unobservable  at 
room  temperature  being  maisked  by  thermal  broadening).  In  general,  any  applied  stress  will 
not  be  imposed  on  the  crystal  along  its  crystallographic  axes,  ais  represented  in  equation  1, 
but  rather  in  some  other  coordinate  frame.  The  stress  on  the  crystal  structure  will  then  be 
given  by  the  applied  stress  components  resolved  onto  the  crystallographic  eoces,  i.e.  through 
an  orthogonal  coordinate  transformation: 

^ij  ~  ^kl  (2) 

where  Cij  is  the  trjinsformation  matrix  relating  the  frame  of  reference  in  which  the  applied 
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stresses  axe  defined  with  respect  to  the  crystallographic  axes  of  the  lattice. 

Thus,  the  frequency  shift  of  a  fluorescence  line  in  a  luminescing  crystal  oriented  at  an 
arbitrary  angle  to  a  superimposed  stress  and  strain  field  is  given  by  the  following  tensoriaJ 
relation: 

—  n,j  aUcOjl  (Tfil  —  n,j  Cijkl  dkmO'ln  ^mn  (3) 

One  of  the  major  limitations  of  piezo-spectroscopy  is  evident  from  this  equation.  Unless 
the  fluorescing  species  emits  at  several  different  frequencies,  so  that  a  set  of  simultaneous 
equations  (one  for  each  transition)  can  be  generated,  there  may  be  an  insufficient  number  of 
distinct  transitions  to  solve  for  each  of  the  maximum  of  six  independent  stress  components 
needed  to  completely  describe  the  stress  state.  In  some  materials,  the  point  symmetry  of  the 
fluorescing  ion  may  also  preclude  the  determination  of  all  the  independent  stress  components. 
This  will  be  seen  to  be  the  case  for  both  MgO  and  AljOa- 

Of  the  many  fluorescence  spectra  reported,  that  of  chromium  fluorescence  in  sapphire 
haa  probably  been  the  most  thoroughly  studied.  The  chromium  ions  replace  aluminum  ions 
substitutionally  in  the  sapphire  lattice  and  maintain  their  trivalent  state,  Cr^'*’.  Although 
the  chromium  ion  is  surrounded  by  an  octahedron  of  oxygen  ions,  the  octahedron  is  distorted 
along  the  body  diagonal  giving  a  point  symmetry  of  D3rf(iJm)  for  the  chromium  ion.  The 
distortion  also  causes  an  asymmetry  of  the  electrostatic  crystal-field  and  results  in  a  splitting 
of  the  energy  levels  by  0.004  eV,  in  turn  giving  rise  to  the  two  distinct  radiative  trcinsition 
lines  in  ruby,  Ri  and  Rj  at  1.790  and  1.794  eV  respectively.  From  the  point  symmetry  of  the 
chromium  ion,  it  can  be  shown  that  the  off-diagonal  components  of  the  piezo-spectroscopic 
tensor  are  zero  and  that  the  diagonad  components  are  symmetrical  about  the  c-aucis,  namely 
rill  =  Hja  1133.  For  simplicity,  we  will  use  On  =  022  =  Ha  and  1133  =  He  for  Ri  and  R2 
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lines.  Substituting  these  components  into  equation  3,  the  frequency  shift  for  a  chromium 
fluorescence  line  can  be  expressed  in  terms  of  the  applied  stress  components  as: 

Al/  =  na((7u  -f  <7ii  +  <Tii)  +  (lie  —  na)(a3i  <Tii  +  a|2<^22  +  alaO'Ss) 

-|-  2(nc  —  na)(a3ia32  <Ti2  +  032^^33  <^23  +  031^33  <^3l)  (4) 

This  differs  from  the  expression  derived  by  Grabner  (equation  4  reference  2)  in  having  an 
additional  term  (the  second)  proportional  to  the  trace  of  the  applied  stress  tensor.  For  the 
particular  case  in  which  the  applied  stress  is  purely  hydrostatic,  this  equation  simplifies  to 
that  used  in  the  high-pressure  physics  and  materials  literature,  namely 

At/  =  -(2na  +  ne)P  (5) 

where  P  is  the  pressure.  Also,  unlike  Grabner’s  equation,  equation  4  reduces  to  equation  1 
when  the  stress  field  is  applied  along  the  crystal  axes.  The  point  symmetry  of  the  chromium 
ion  simplifies  equation  3  but  nevertheless  the  stress  components  cannot  be  determined  com¬ 
pletely  from  the  shifts  of  the  Ri  and  R2  lines  unless  the  misorientation  between  the  applied 
stress  and  the  crystal  axes  is  known  from  some  other,  independent  meeisurement. 

The  chromium  fluorescence  spectrum  from  magnesia  has  also  been  investigated  in  some 
detail  [7],  and  can  similarly  be  used  to  monitor  strains.  As  in  sapphire,  the  chromium  ion 
occupies  substitutionally  a  lattice  site  in  magnesia  and  is  surrounded  by  an  octaJiedron  of 
oxygen  ions.  However,  the  oxygen  octahedra  remain  undistorted  so  the  chromium  point 
symmetry  is  raised  to  m3m.  As  a  consequence,  the  level  is  not  split  and  only  one  line, 
at  14319  cm“\  occurs.  The  m3m  point  symmetry  imposes  the  additional  requirement  that 
the  diagonal  piezo-spectroscopic  coefficients  are  all  equal  (Hu  =  n22  =  Hsa).  Thus,  for 
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chromium  luminescence  in  magnesium  oxide,  the  line  shift  is  sensitive  only  to  the  trace  of 
the  stress  tensor: 

Ai/  =  +  (7u  -I- <7ii)  (6) 

This  insensitivity  to  the  deviatoric  stresses,  nevertheless,  has  the  advantage  that  the  shift  is 
independent  of  the  orientation  of  the  crystal  structure  with  respect  to  the  applied  stress. 

Rare  earth  ions,  such  as  Sm^"*"  and  Nd^"*",  can  substitute  for  the  yttrium  ion  in  the  yttrium 
aluminum  oxide,  garnet  structure  and  give  rise  to  sharp  characteristic  fluorescence  spectra. 
The  three  possible  substitutional  sites,  although  equivalent  under  the  cubic  point  group 
symmetry  of  YAG,  have  orthorhombic  site  symmetry.  Since  ions  at  all  three  sites  will 
fluoresce  when  the  crystal  is  excited,  applying  a  stress  should  in  general  breedc  the  degeneracy 
associated  with  the  equivalence  of  the  three  sites  and  cause  line  splitting.  However,  in 
practice,  the  stress  splitting  is  masked  by  thermal  broadening  at  room  temperature,  and 
so  the  fluorescence  appears  to  come  from  a  crystal  having  cubic  symmetry  [8]  with  equal 
piezo-spectroscopic  coefficients. 

3  INHOMOGENEOUS  STRAIN  FIELDS 

The  above  analysis  pertauns  to  the  case  of  a  single  crystal  subject  to  a  homogeneous  strain 
field  over  the  volume  of  material  exposed  to  the  excitation  probe.  When  the  probed  volume 
contains  more  than  one  crystallite  or  the  strains  vaiy  within  it,  the  observed  spectral  shifts 
are  more  complex.  They  can,  however,  be  understood  from  simple  superposition  arguments 
provided  inelastic  scattering  of  the  luminescence  is  not  appreciable  and  the  strains  do  not  vary 
over  distances  comparable  to  the  lattice  spacing.  (For  the  case  of  the  N-  and  F-  fluorescence 
lines,  the  latter  requirement  is  a  little  more  stringent.  In  that  case,  the  strains  should  not 
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vary  over  distances  associated  with  the  separation  of  the  ions  making  up  the  chromium  pairs 
from  which  these  characteristic  lines  originate). 

Subject  to  these  restrictions,  strain  inhomogeneities  and  the  presence  of  multiple  grains 
within  the  probed  volume  can  cause  both  broadening  and  shifting  of  the  luminescence  lines 
as  analyzed  in  Appendix  A.  In  essence,  each  element  of  material  over  which  there  is  a 
constant  crystedlographic  orientation  and  homogeneous  streiin  contributes  a  signal  intensity 
proportional  to  its  volume  and  concentration  of  fluorescence  ions.  The  spectral  shift  is  given 
by  equation  3  with  the  appropriate  values  for  each  of  the  parameters.  The  superposition  of 
the  lines,  shifted  according  to  the  strain  of  each  volume  element,  causes  the  observed  line 
to  appear  to  be  broadened  and  shifted.  The  broadening  is  not  a  measure  of  the  deviatoric 
strains  per-se  as  has  been  stated  in  the  literature  previously  [2]  but  rather  of  the  statistical 
variation  in  strains  and  crystallographic  orientations  in  the  region  examined.  As  a  result  the 
broadening  need  not  be  symmetrical.  Indeed  the  broadened  lines  can  have  a  fine  structure 
resulting  from  there  being  one  or  two  regions  whose  contributions  are  disproportionately 
l2irge.  This  is  probably  the  cause  of  the  fine  structure  in  the  spectra  reported  by  Grabner 
[2]. 

Although  analytical  relationships  between  line  shifts,  broadening  and  inhomogeneous 
stress  state  in  polycrystalline  materials  can  not,  in  general,  be  presented,  a  number  of  simple 
cases  of  practical  importamce  can  be  considered.  Equations  for  a  number  of  these  are  derived 
as  described  in  the  following  sub-sections. 

3.1  Polycrystalline  Material 

The  first  example  for  which  an  analytical  function  can  be  derived  corresponds  to  a  poly- 
crystadline  material,  with  no  crystallographic  texture  or  internad  residual  stress,  under  a 
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homogeneous  applied  stress  and  probed  with  an  excitation  beam  large  compared  to  the 
grain  size.  In  this  case,  the  contributions  from  the  individual  crystallites  are  averaged  over 
many,  randomly  oriented  grains. 

The  frequency  shift  in  the  fluorescence  line  is  given  by  the  first  moment  of  the  distribution 
of  fluorescence  frequencies.  Using  equations  3,  and  A6  and  B3  in  the  appendices,  this  can 
be  simply  expressed  as: 

Ai/  =  J J J  P{6,<f>,il;)Avd6d<i>dil; 

which  reduces  to 

Au  =  —  (Till  -h  1122  +  n33)(<Tii  +  <Tii  +  Cii)  (7) 

This  equation  shows  that,  irrespective  of  the  crystal  structure,  the  frequency  shift  of  the 
fluorescence  from  a  large  number  of  grains  in  a  polycrystalline  material  is  dependent  only 
on  the  trace  of  the  applied  stress  tensor. 

The  polycrystalline  ensemble  of  grains  also  leads  to  an  intrinsic,  stress  dependent  broad¬ 
ening  of  the  fluorescence  line.  The  peak  broadening  is  given  by  the  second  moment  of  the 
distribution,  which  using  equation  A8,  may  be  written  ats: 

<  Av^  >  =  —  |(n22  —  +  (Has  —  IIn)^  —  ((1122  —  nii)(n33  —  Hu)] 

X  -  (o'?!  4-  (722  ^33  ~  —  <^22<^33  —  <733<7ll)  +  (cTjj  -j-  <733  +  <7|i)j  (8) 

Such  broadening  is  due  to  the  differences  in  the  diagonal  piezo-spectroscopic  coefficients. 
Consequently,  for  matericils  with  fluorescing  ions  having  cubic  point  symmetry,  such  ais  Cr^"*" 
in  MgO,  (for  which  the  piezo-spectroscopic  coefficients  are  equal)  the  application  of  an  ex¬ 
ternal  stress  does  not  cause  any  broadening. 


8 


3.2  Inhomogeneously  Stressed  Single  Crystal 

The  second  case  corresponds  to  a  single  crystal  containing  inhomogeneously  stressed  re¬ 
gions.  A  special  case,  but  nevertheless  am  important  one,  is  one  in  which  the  stress  varies 
with  depth  below  the  surface,  such  as  might  pertain  to  surface  stressing  produced  by  ion 
exchange  or  when  the  sample  is  loaded  in  bending.  For  purposes  of  generality,  we  may  as¬ 
sume  that  the  concentration  of  the  fluorescing  ion,  c,  is  a  function  of  depth  z.  We  can,  in 
principle,  also  incorporate  various  instrumental  functions,  such  as  the  collection  efficiency 
as  a  function  of  depth  of  the  apparatus  used  to  measure  the  fluorescence  spectrum,  into  a 
function  g{z).  The  function  g{z)  allows  the  fluorescence  signal  to  be  normalized: 

J  c{z)g{z)dz  =  1 

and  a  weighting  function  introduced: 

W{^u)di/  =  c{z)g{z)dz  (9) 

Again,  for  purposes  of  generality,  the  frequency  shift  can  be  written  so  «is  to  include  a  shift 
dependent  on  the  concentration,  c,  of  the  fluorescing  species: 

Ai/ =  Ili,a,kajiaki  + I3{c- Co)  (10) 

where  cq  is  the  concentration  of  the  reference  sample  which  the  shift  is  compared. 

The  average  peak  shift,  expressed  by  the  first  moment  of  the  fluorescence  peak  then 
becomes 

^=/[n  ij  OikOji  (Tki  +  I3{c  -  Co)]  c{z)g{z)dz  (11) 
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The  broadening  is  a  mixture  of  the  inhomogeneous  broadening  due  to  the  depth  depen¬ 
dence  of  shift  and  the  intrinsic  homogeneous  broadening  due  to  change  of  concentration.  For 
simplicity,  we  give  the  result  for  the  case  of  constant  concentration  cq: 

<  Aj/^  >=  J {Uijaikajiakiyc{2)g{z)dz  -  (12) 

in  which  ST  ^  is  also  evaluated  at  c  =  cq. 

3.3  Anisotropic  Thermal  Expansion  Residual  Stresses 

The  third  case  considered  here  is  of  a  polycrystalline  materiaJ  for  which  the  thermal  ex¬ 
pansion  coefficient  is  crystallographic<illy  anisotropic.  As  a  result  of  such  thermal  expansion 
anisotropy,  residual  stresses  can  be  created  on  cooling  from  the  fabrication  temperature,  with 
the  magnitude  of  the  stresses  being  dependent  on  the  geometrical  constraint  imposed  by  the 
neighboring  grains,  the  elastic  properties  of  the  grains  emd  the  extent  to  which  accommo¬ 
dation  mechanisms,  such  as  creep  and  plastic  deformation,  which  cam  relax  the  stresses.  In 
the  special  cases  in  which  the  grains  do  not  constrain  one  another  on  cooling  or  the  accom¬ 
modation  mechanisms  are  complete,  no  internal  stresses  develop  and  neither  a  line  shift  nor 
broadening  would  arise.  However,  there  is  in  genertil  some  mutucil  constreunt  imposed  by  the 
contracting  grains  that  causes  stresses  to  develop  within  the  material.  The  average  stress 
over  the  material  must  be  zero  but  variations  in  stress  from  one  grain  to  another  can  cause 
both  a  broadening  of  the  line  and  a  line  shift.  Although  the  specific  values  will  vary  from 
place  to  place  in  the  sample,  since  the  stresses  will  also  vary  from  grain  to  grain,  it  is  possible 
to  calculate  the  average  values  as  will  now  be  shown.  Two  assumptions  are  mawie.  First, 
that  the  probed  volume  encompasses  a  sufficiently  large  number  of  grains  that  an  ensemble 
average  can  be  considered,  and  secondly  that  the  distribution  of  stresses  can  be  described 
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by  an  appropriate  distribution  function. 

Consider  a  material,  such  as  alumina,  which  has  a  thermal  expansion  coefficient  that  is 
larger  along  its  c-axis  than  along  its  a-axis,  and  has  an  axis  of  rotational  symmetry  along 
the  c-axis.  On  cooling  from  its  fabrication  temperature,  tensile  stresses  will,  on  average, 
develop  aJong  the  c-axis  of  the  grains  and  compressional  stresses  will  develop  in  the  bas2d 
plane.  Therefore, 

37  <  0  ,  57  >  0  (13) 


Now,  consider  a  grain  intersected  by  an  arbitriurily  placed  plane  that  divides  the  poly¬ 
crystalline  sample  into  two.  The  orientation  of  the  graun  can  be  specified  by  the  angular 
coordinates  (^,  <^)  of  its  c-axis  (figure  1).  Provided  the  dividing  plane  intersects  a  large  num¬ 
ber  of  grains,  with  all  possible  stress  states,  we  need  only  consider  the  average  stresses  along 
the  a-  and  c-  axes.  Therefore  the  stress  in  the  basal  plane  is  equi-biaxial,  and  rotation  of 
the  grain  around  the  c-axis  does  not  alter  its  stress  state.  The  force  acting  per  unit  area  in 
any  grain  perpendicular  to  the  dividing  plane  can  be  expressed  by  resolving  the  stresses  in 
that  grain: 

/  =  ^sin^ -I- o^cos^  (14) 


Since,  across  the  dividing  plane,  the  orientation  of  the  grains  is  random,  the  average  force 
Ccin  be  written  as: 

-  1  /■’f/2 

f  =  7r  /sin^de  (15) 

Jo  Jo 


Since  the  force  must  on  average  be  zero,  /  =  0,  the  mean  stresses  in  the  c-  and  a-  crystallo¬ 
graphic  directions  are  related  by: 


<Tc  = 


(16) 
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Recognize  that  the  average  stresses  described  above  are  equivadent  to  volume  average,  the 
mean  shift  of  the  fluorescence  line  is  therefore, 


Ai/  =  ^  Ha^  (Tc 


The  peak  broadening  can  be  calculated  using  equation  A8  on  the  assumption  that  the 
distribution  in  stresses  in  the  a-  and  c-  crystallographic  directions  are  independent  of  one 


auiother: 


P{<Ta,<rc)  =pMp{<Tc) 


and  are  also  Gaussian  with  the  same  width,  <7^ 


.  .  1  (<Ta-<Ta)2 

2ai 


The  assumption  of  independence  of  the  stress  distributions  in  the  different  directions  enables 
the  weighting  function  in  equation  A8  to  be  written  as: 


W{Al/)dv  =  P{<Ta,<Tc)d<Tad(Tc  =  p{a  a)p{cr  c)d(T  ^dc  c 


So,  the  peak  broadening  becomes: 


<  Ai/*  >=  411^  (tI  +  n*  cr^  +  4nanc  CTc 


Thus,  using  measurements  of  the  mean  peak  shift,  which  from  equation  17  gives  the  mean 
stresses  along  the  a-  and  c-  axes,  and  a  measurement  of  the  peak  width,  the  stress  distribution 
within  the  material  can  be  determined. 
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These  equations  shows  that  the  crystallographic  anisotropy  in  thermal  expansion  leads 
to  both  am  average  shift  and  broadening  of  the  fluorescence  line.  They  are  related  to  the 
average  values  of  the  stresses  and  the  stress  distribution.  The  spatially  varying  stresses  can 
be  significantly  larger  than  the  average,  especially  those  acting  over  smaller  distances  than 
the  grain  size,  related  to  the  presence  of  the  grain  edges  and  corners.  These  are  primarily 
responsible,  for  instance,  for  micro-crack  extension.  Although  detailed  information  of  such 
stresses  does  not  appear  to  be  amenable  to  the  analysis  of  the  shape  of  the  fluorescence  lines, 
their  magnitude  amd  the  occurring  frequency  are  reflected  by  the  stress  distribution. 

3.4  Matrix  Stresses  Due  to  Presence  of  Fibrous  Inclusions 

The  final  case  considered  is  that  of  a  polycrystalline  matrix  surrounding  an  embedded 
fiber.  The  radial  and  tangential  stresses  in  a  matrix  around  a  misfitting  fiber  of  radius,  tq, 
are  given  by  standard  elastic  solutions,  varying  with  the  inverse  square  of  distance,  r,  from 
the  fiber  axis.  For  observations  made  in  a  cross-section  perpendicular  to  the  fiber  axis  (axis 
3),  the  matrix  stresses  may  be  written  as 

<Tii  =  —<T22  =  ;  <733  =  0  (21) 

Thus,  from  equation  3,  the  frequency  shifts  are  given  by 

A»/  =  (1122  ~  nu)(a2i  —  Ojj)  (Til  -I-  (1133  —  nii)(a3i  —  032)  <T22  (22) 

which  would  be  detectable  with  a  probe  having  high  optical  resolution  (<C  tq).  If  the  probe 
is  large,  such  that  either  many  grains  in  a  polycrystalline  matrix  or  regions  within  a  circle 
encompassing  the  fiber  in  a  single  crystal  matrix  aie  excited  simultaneously,  then  the  average 
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peak  shift  is  Ai/  =  0.  The  peak  broadening,  for  such  a  large  probe,  would  also  average  out 
for  single  crystal  and  polycrystalline  material  to  a  v£ilue  of 

<  Ai/^  >=  —  <7ii  [(1122  —  Hii)^  +  (Haa  —  —  (1123  ~  nu)(na3  —  Hn)]  (23) 

where  <7^^  is  averaged  over  the  volume  probed. 

4  OBSERVATIONAL  TESTS 

In  this  section  a  series  of  measurements  are  described  that  provide  experimental  tests  and 
validations  of  the  equations  presented  in  the  preceding  sections.  The  tests  were  made  with 
aluminum  oxide  because  of  its  ready  availability.  The  experiments  were  performed  using  an 
unmodified  optical  microprobe‘  to  both  excite  the  luminescence  ^lnd  to  collect  and  analyze 
the  resulting  fluorescence  spectrum  using  an  attached  spectrometer.  (The  optics  of  the 
Microprobe  have  been  described  in  an  ewlier  publication  [9]).  A  number  of  laser  frequencies 
have  been  used  to  excite  the  luminescence  but  the  results  reported  here  were  all  obtained 
using  an  argon  ion  laser  with  a  wavelength  of  514  nm.  The  experimental  procedure  was  to 
select  a  region  of  interest  on  the  sample  with  the  attached  optical  microscope,  focus  the  laser 
beam  to  a  spot  on  the  selected  feature  and  analyze  the  fluorescence  produced.  The  intensity 
of  the  Ri  and  R2  fluorescence  lines  were  typically  scanned  by  integrating  over  0.5  second 
intervals  at  a  spacing  of  0.2  wavenumbers,  with  the  intensities  being  recorded  under  computer 
control.  The  collected  data  was  subsequently  analyzed  with  curve  fitting  algorithms,  included 
in  the  LabCalc  software  package^,  to  identify  the  position  of  the  fluorescence  peak.  By  using 
objective  lenses  of  50  x  and  100  x  magnifying  powers  minimum  spot  sizes  of  ~  5^m  and 

^  InatrumcnU  SA,  Raman  Microproba  Model  T64000. 

^Gatictk  Indiiatriea  Corp. 
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1  —  2fim  can  be  produced.  Larger  probe  sizes  were  formed  by  partial  defocusing.  Also, 
by  appropriate  choice  of  the  collection  aperture  size,  the  attainable  axial  resolution  could  be 
varied  by  the  collection  optics  from  lO/xm  to  ~  200/im  [10]. 

All  the  measurements  were  made  with  the  sample  at  room  temperature.  Although  a 
change  in  temperature  is  known  to  cause  a  line  shift  of  ~  0.14  cm~^/°C  [11]  (to  smaller 
wavenumbers  with  increasing  temperature),  no  line  shift  was  noted  when  the  spot  size  was 
systematically  decreased  (thereby  increasing  the  power  density  incident  on  the  sample)  in¬ 
dicating  that  no  significant  heating  of  the  probed  volume  occurred.  This  lack  of  sensitivity 
to  spot  size  was  attributed  to  both  the  low  absorption  of  ruby  at  the  excitation  frequencies 
and  the  relatively  low  illumination  intensities  used.  Variations  were,  however,  noted  with 
variations  in  room  temperature  and  so  temperature  corrections  were  made  to  the  peak  shifts. 

4.1  Polycrystalline  Sample  Under  Stress 

In  order  to  test  the  validity  of  equations  7  and  8  concerning  the  stress  dependence  of  the 
peak  shift  and  broadening  in  a  polycrystalline  material,  polished  bars  of  polycrystalline  alu¬ 
mina  (Coors  AD96)  were  subject  to  an  applied  stress  and  the  fluorescence  spectra  recorded. 
The  stress  was  applied  using  a  four-point  bend  stage  specially  designed  to  fit  under  the  op¬ 
tical  microscope  of  the  microprobe.  The  polycrystalline  alumina  contained  sufficient  levels 
of  chromium  impurity  to  provide  sufficiently  strong  fluorescence  peaks  to  be  analyzed.  In 
order  to  average  over  many  grains,  a  probe  diameter  of  ~  50/im  was  used  for  the  fluorescence 
measurements.  The  grain  size  of  the  alumina  was  2-20/im. 

Since  alumina  is  thermally  2inisotropic,  the  analysis  in  section  III.l  has  to  be  extended  to 
separate  the  effect  of  thermal  residual  stress  existing  in  the  sample  from  any  applied  stress. 
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The  shift  from  a  small  volume  element  dv  (much  smaller  than  the  grain  size)  is 


Ai/  =  2na(<T^  +  +  nc(<T^  +  (tJ  )  (24) 

where  the  superscripts  A  and  T  denote  applied  «md  thermal  stresses  respectively.  Since  the 
above  equation  is  linear,  the  effects  of  two  types  of  stress  can  be  separated; 

Ai/  =  Ai/^  +  Aj/^  (25) 


Therefore,  the  mean  shift: 

(  Aj/^  +  Au^)dv  =  +  A^  (26) 

Since  the  shift  generated  by  thermal  stress  (see  section  4.3)  can  be  measured  independently, 
Ai/^  can  be  obtained. 

The  broadening  is  calculated  by  using  equation  A8: 

<  Au^  >  =  ^  +  Ai/'^fdv -Au^ 

=  <  {Av^)^  >  +  <  {Au'^f  >  +2  [Au'^AuT  -  Ai^  AJ^]  (27) 

So  the  broadening  can  be  separated  provided  that  Ai/‘^  and  Ai/^  are  uncorrelated. 

The  measured  peak  shifts  for  the  Ri  and  R2  fluorescence  lines  as  a  function  of  applied 
stress  are  shown  in  figure  2.  The  lines  through  the  data,  fitted  by  a  least  square  analysis, 
shows  that  the  stress  dependence  of  the  shift  of  the  Ri  and  Ra  fluorescence  lines  is  2.46  and 
2.50  cm“*GPa~*  respectively.  The  fit  coefficient  R  is  0.997  for  both  lines.  For  comparison 
to  the  theory  in  section  III.l,  equation  7  would  predict  that  when  a  uniform  stress  is  applied 
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to  polycrystalline  alumina,  the  peak  shift  would  be  given  by: 


Ai/  =  i(2n,  +  n,)<T^ 

Using  the  hydrostatic  pressure  dependence  of  peak  shifts  for  the  Rj  and  R2  lines  given  by 
Munro  aL  [12],  the  sum  of  the  spectroscopic  coefficients,  2no  +  lie  have  the  values  of 
7.59  and  7.62  cm~^  GPa~^  for  the  Rj  and  Rj  lines  respectively.  The  peak  shifts  of  the 
Ri  and  Rj  lines  would  then  be  2.53  and  2.54  cm“^GPa~*,  close  to  that  observed.  The 
uniaxial  compression  experiment  by  Feher  and  Sturge  [13]  provided  values  for  Ho  and  He 
independently  with  values  of  2.7  and  2.15  cm~*  GPa“*  for  both  Ri  and  R2  lines.  Using 
this  data,  the  predicted  stress  dependence  would  be  2.52  em'^GPa"',  again  close  to  that 
observed.  We  can  therefore  conclude  that  our  measured  shifts  fall  within  the  uncertainty  of 
the  existing  data  on  the  piezo-spectroscopic  coefficients. 

The  measured  peak  widths  for  the  Ri  and  R2  fluorescence  lines  as  functions  of  applied 
stress  are  shown  in  figure  3.  For  this  material  the  systematic  broadening  due  to  the  applied 
stress  of  both  lines  is  small  compared  to  the  variability  from  one  point  on  the  sample  to 
another. 

4.2  Single  Crystal  With  Surface  Stress 

A  single  crystal  treated  to  produce  a  surface  stress  represents  a  special  case  of  the  inhomo- 
geneously  stressed  single  crystal  considered  in  section  III. 2.  To  provide  an  experimental  test 
of  the  effect  of  inhomogeneous  stressing,  a  residuail  surface  stress  was  produced  in  a  single 
crystal  of  c-axis  sapphire  by  in-diffusion  of  chromium  and  the  fluorescence  meaisured  as  a 
function  of  an  applied  stress.  Since  chromium  oxide,  Cr203,  is  isostructural  with  AI2O3  but 
has  a  larger  lattice  parameter,  diffusion  of  chromium  into  the  surface  of  AI2O3  is  expected 
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to  produce  a  compressive  surface  stress. 

As  chromium  is  both  the  fluorescing  ion  and  the  ion  causing  the  surface  stress,  the  analysis 
of  section  III. 2  must  include  the  stress  created  by  the  in-diffusion  of  chromium.  In  this  case, 
the  frequency  shift  is  expected  to  be  the  sum  of  three  distinct  contributions,  that  due  to 
«uay  applied  stress,  that  due  to  diffusion-induced  stresses  <7^  and  that  due  to  variations  in 
fluorescing  ion  concentrations: 


Au  =  Ilij  <7,^  -I-  Ilij  (7,5  +  0{c  -  Co)  (28) 

The  overall  pe2J{  shift  from  the  in-diffused  layer  can  be  written  as: 

Ai/  =  J  An  c(z)g{z)dz  =  Ai^^  +  Av^  -{■  Av^  (29) 

Since  the  frequency  shift  in  equation  28  is  the  sum  of  the  individual  components,  the  mean 
peak  shifts  of  each  component  are  also  linearly  independent  as  shown  by  equation  29.  Thus, 
in  zm  experiment  in  which  am  external  stress  is  applied  to  a  residually  stressed  surface,  the 
mean  frequency  shift  observed  is  predicted  to  be  due  to  the  applied  stress  only: 

Ai/^  =  j  Av^ c{z)g{z)dz  =  11, ^  cr,^  (30) 

For  four-point  bending,  the  applied  stress  is  essentiadly  constant  over  the  depth  of  the 
in-diffused  layer,  and  so  the  mean  frequency  shift  will  be  linearly  related  to  the  applied 
stress. 

For  the  experiment,  in-diffusion  of  chromium  was  performed  in  a  three  step  process.  In 
the  first  step,  a  chromium  acetate  solution  was  spun-coated  onto  a  c-eixis  oriented  sapphire 
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disc.  The  coating  was  then  pyrolyzed  and  oxidized  by  heating  to  600°  C  in  aur.  Then, 
finally,  the  disc  was  annealed  at  1600°C  to  diffuse  in  the  chromium  and  beams  cut  from  the 
disc.  To  investigate  the  effect  of  applied  stress  the  sapphire  beam  was  loaded  in  four-point 
bending  with  c-axis  perpendicular  to  the  loading  plane,  placed  under  the  microscope  of  the 
microprobe,  and  the  resulting  fluorescence  recorded  as  a  function  of  load.  The  observed  peak 
shift  for  the  Ri  and  R2  lines  as  a  function  of  applied  stress  is  shown  in  figure  4.  For  this 
data,  the  stress  dependence  is  2.42  and  2.61  cm~^GPa”^  for  the  Ri  and  Rj  lines  respectively, 
which  compares  favorably  to  the  mean  Ila  value  of  2.7  cm“^GPa“^  given  in  [13].  The  fit 
coefficient  R  equals  0.996  and  0.994,  respectively. 

4.3  Thermal  Expansion  Anisotropy  Stresses 

Residual  stresses  are  known  to  be  produced  in  polycrystalline  alumina  on  cooling  as  a 
result  of  the  crystallographic  anisotropy  in  thermal  expansion  coefficients  along  the  c-  and 
a-  axes.  In  order  to  measure  the  average  residual  stresses,  the  fluorescence  from  three 
different  sintered,  polycrystalline  aluminas  (Coors  AD96,  Coors  AD  995  and  UCSBl)  and  a 
polycrystalline  alumina  hipped  at  UCS3  (UCSB2)  was  recorded  using  an  excitation  probe 
diameter  of  50/im.  The  UCSBl  material  is  a  high  purity  alumina  sintered  at  1800°C  for  1  hr. 
in  hydrogen.  The  probe  size  was  much  larger  than  the  grain  size  of  the  ceramics,  assuring 
that  the  fluorescence  was  being  averaged  over  a  large  number  of  grains.  The  ceramic  samples 
contained  an  adequate  chromium  impurity  level  for  the  Ri  and  R2  fluorescence  pealcs  to  be 
recorded  with  sufficient  signal- to- background  ratio  that  precise  measurements  could  be  made 
of  the  peak  shape  and  peak  shift.  The  peak  shift  and  peaJc  broadening  were  measured  relative 
to  those  recorded  from  a  polished  single  crystsJ  of  sapphire  that  also  contained  a  trace  level 
of  chromium  impurity.  For  comparison,  the  spectrum  obtained  from  the  AD995  sample  is 
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plotted  together  with  that  of  a  single  cryst£il  sapphire  in  Figure  5.  Several  spectra  were  taken 
at  different  spots  for  each  sample,  the  average  peak  shifts  and  broadening  for  the  samples 
axe  listed  in  table  I.  Also  shown  in  the  table  axe  the  average  stresses,  and  the  breath 

of  stress  distribution  cr^  calculated  using  equations  16,  17  and  20  above. 

5  CLOSING  REMARKS 

The  preceding  analyses  and  experiments  illustrate  how  the  technique  of  fluorescence  spec¬ 
troscopy  can  be  used  to  make  measurements  of  commonly  arising  types  of  stresses  in  ceram¬ 
ics.  This  is  possible,  despite  the  fact  that  the  fluorescence  shift  is  a  scalar  quantity  whereas 
the  complete  stress  state  is  tensorial,  because  considerable  simpliflcation  of  the  equations 
(equation  3)  can  be  made  in  many  practical  situations.  For  instance,  stresses  can  be  applied 
at  prescribed  orientations  to  the  crystallographic  axes  of  single  crystals.  Also,  the  poly- 
crystaJline  nature  of  many  ceramics  enables  the  effects  of  crystallographic  orientation  to  be 
averaged  out.  For  instance,  as  equation  7  demonstrates,  the  frequency  shift  from  a  stressed 
polycrystalline  ceramic  is  independent  of  its  crystal  structure  and  is  only  dependent  on  the 
trace  of  the  applied  stress. 

Two  of  the  analytical  results  are  of  particular  interest.  One,  is  the  finding  that  the 
fluorescence  shift  is  dependent  only  on  the  applied  stress  and  not  on  the  residual  stresses 
produced  by  either  thermal  expainsion  mismatch  «is  in  the  zase  of  polycrystailine  AI2O3,  or 
by  ion  in-diffusion  as  in  the  case  of  Cr^'*'  in-diffusion  into  sapphire  single  crystal  .  The  later 
also  suggests  that  a  fluorescent  dopant  ion  can  be  deliberately  introduced  into  a  structure  or 
materiid  to  monitor  the  applied  stress  through  changes  in  the  fluorescence  frequency.  In  fact, 
we  have  been  utilizing  this  cheiracteristics  to  make  measurements  of  stresses  in  metal-ceramic 
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laminates  [14]. 

The  second  notable  result  concerns  the  me«isurement  of  the  residual  stresses  due  to 
anisotropic  thermal  expansion.  The  analysis  indicates  that  the  mean  stress  in  the  a-  and 
c-  crystallographic  directions  in  aluminum  oxide  can  be  obtained  from  the  mean  frequency 
shift,  cind  that  the  distribution  in  internal  stresses  can  also  be  determined.  These  results  rep¬ 
resent  an  advamce  on  previous  studies  of  anisotropic  thermal  expansion  stresses  in  aluminum 
oxide  [2,15].  As  already  mentioned,  Grabner’s  work  led  to  am  incomplete  expression  for  the 
relationship  between  the  frequency  shift  and  stress,  and  so  the  values  he  derived  axe  probably 
unreliable.  The  work  of  Blendell  and  Coble  [15],  was  primarily  concerned  with  the  dilFusional 
relaxation  of  the  anisotropic  thermal  expansion  stresses.  They  assumed  an  average  stress  in 
the  material  and  calculated  its  relaxation  for  different  grain  sizes  and  cooling  rates.  They 
did  not  relate  this  average  stress  to  any  of  the  crystallographic  stresses  nor  did  they  describe 
how  they  determined  this  average  stress  from  their  fluorescence  observations.  However,  the 
measurements  in  table  I  indicate  that  the  mean  stresses  in  the  samples  of  polycrystalline 
alumina  examined  are  similar  in  magnitude  to  those  calculated. 
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APPENDIX 


A  Peak  Shifts  and  Broadening  due  to  Stresses 

If  the  probed  volume  were  strain  free,  the  fluorescence  peaks  would  be  Lorentzian  in  shape 
and  have  a  breadth  dependent  only  on  the  finite  temperature  of  the  sample.  The  presence 
of  stresses,  however,  causes  the  peaks  to  be  shifted  and  dispersed  as  shown  in  the  following. 

Let  the  peak  shape  in  the  unstrained  crystal  be  represented  by  the  frequency  function, 
/(»/),  with  the  peak  centroid  given  by: 

^0  =  j  f{v)vdv  (Al) 

and  the  area  of  the  peak  is  normalized  to  be  unity: 

j  =  1  (A2) 

The  dispersion  of  the  peak,  corresponding  to  the  peak  width,  can  be  expressed  by  the 
second  moment  of  the  distribution: 

<  >0=  j  f{v)[y  -  (A3) 

The  elfect  of  amy  stress  distribution  can  be  related  to  meaisured  peak  shape  by  assuming 
that  the  probability  of  causing  a  frequency  shift,  At',  can  be  expressed  as  a  weighting  function 
W{^v).  The  resulting  peak  is  then  given  by  the  function; 

F{v)  =  J  fii'-  ^u)W{Au)dAu  (A4) 
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Since  /  W{Av)dAi/  =  1,  /  F{u)dv  =  1  the  axea  under  the  pe£ik  is  conserved. 

The  centroid,  or  mean  value,  of  the  peak  from  the  stressed  volume  is  given  by  the  first 
moment  of  the  function  representing  the  peak: 

=  J  F{v)i>dv  (A5) 

The  peak  shift  is  then: 

Av  =  Vc  —  I'o 

=  J  Ai^W{Au)dAu  (A6) 

As  spectra  /(t/)  and  F{i/)  can  be  obtained  from  unstressed  and  stressed  regions,  the  mean 
shift  Ai/  can  be  determined  and  from  that  the  stress  deconvoluted. 

The  width  of  the  stress  broadened  peak,  F(i/),  is: 

<  1/2  >p=  j  F{i>){iy  -  Vc)^di/  (A7) 

The  broadening  is  then  seen  to  be  the  difference  between  the  measured  widths  of  the  peaks 
from  the  strain  free  regions  and  from  the  stressed  regions: 

<  1/2  >f  -  <  1/2  >0  =  j  Av^  W{Ai')dAi/  —  Av^ 

=  <  >  (A8) 

Again,  the  measured  broadening  and  the  stress  distribution  is  related  by  the  weighting 
function. 
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B  Polycrystalline  Averages  of  Direction  Cosines 


Collected  here  for  convenience  are  the  polycrystalline  averages  of  the  direction  cosines  a,_, 
appearing  in  tensor  transformation  relations  of  the  form: 

^ij  ~  O’imOjn  (^mn  (®1) 

In  terms  of  the  usuaJ  Euler  angles  {6,  <^,  ip),  the  transformation  matrix  is 

cos  cos  —  sin  cos  0  sin  —  cos  (p  sin  ip  —  sin  (j>  cos  9  cos  ip  sin  0  sin  ^ 

sin  cos  V’ +  cos  cos  sin  —  sin  sin  0  +  cos  cos  ^  cos  0  —  sin^cosc^ 

sin  9  sin  ip  sin  9  cos  ip  cos  9 

where 

0<<p<2ir,  0<9<ir  ,  0<ip  <2-k  (B2) 

If,  as  in  a  polycrystalline  material,  the  orientations  of  the  grains  are  at  random,  the 
probability  a  grain  director  is  oriented  at  {9,  <p,  ip)  is  given  by 

P{9,<p,ip)  =  sin  9d9d<pdip  (B3) 

OTT^ 

By  averaging  over  all  possible  angles,  the  following  relations  used  to  derive  equation  7 
and  8  are  obtained: 


^ij^ik  —  0  1  7  ^  ^  i 
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-TT  2 
“0°/*  =  T7  ’ 


j  ; 


1 


=  l5  . 


j  ] 


^21^22023  —  a2ia22<*23  =  021022^23  =  ^  ; 


and 


a2l022<J31<232  —  <*21<*23fl3lfl33  =  022023032^33  = 


J_ 
15  ■ 
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FIGURE  CAPTIONS 


1.  Schematic  diagram  of  a  representative  grain  cut  by  an  arbitrary  dividing  plane  (^  =  0) 
through  a  polycryst2Llline  material,  illustrating  the  coordinates  and  angles  used  in  the  text 
in  describing  the  orientation  of  the  grain. 

2.  Measured  peak  shifts  of  the  Ri  (top)  and  Rj  (bottom)  fluorescence  lines  from  a 
polycrystedline  alumina  sample  as  a  function  of  applied  stress.  The  sample  was  loaded 
under  four-point  bending  and  spectra  were  recorded  from  different  positions  on  both  the 
tensile  and  compressional  sides  and  over  the  stress  range  indicated.  The  solid  line  through 
the  data  points  represents  the  line  of  best-fit.  Its  slope  is,  within  experimental  uncertainty, 
the  same  as  that  predicted  using  equation  7. 

3.  Measured  peak  broadening  of  the  Rj  and  Rj  fluorescence  lines  as  a  function  of  stress. 
The  broadening  was  obtained  from  the  same  spectra  as  were  used  to  make  the 
measurements  in  figure  2. 

4.  Peak  shift  of  the  Ri  and  Rj  fluorescence  lines  as  a  function  of  applied  stress  for  a 
sapphire  crystal  containing  a  surface  stressed  layer.  The  solid  line  is  the  line  of  best  fit  to 
the  data.The  surface  stress  was  produced  by  in-diffusion  of  chromium.  Details  are  given  in 
the  text. 

5.  A  comparison  of  the  Ri  and  Rj  fluorescence  lines  from  sapphire  and  a  polycrystaJline 
alumina  (Coors  AD995).  The  broadening  of  the  lines  from  the  polycrystalline  alumina  is  a 
result  of  the  distribution  of  residual  thermal  expansion  anisotropy  stresses  within  the 
material.  Both  materials  contain  a  trace  of  chromium  impurity,  sufficient  to  cause 
fluorescence. 
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TABLE  I 


SAMPLE 

Oc 

yj <  > 

cm“^ 

MPa 

MPa 

cm~^ 

MPa 

AD96 

-0.21 

-100 

160 

2.5 

440 

AD995 

-0.15 

-74 

120 

3.6 

620 

UCSBl 

-0.15 

-74 

120 

1.9 

340 

UCSB2 

-0.23 

-110 

180 

1.4 

270 
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Acoustic  emission  associated  with  the  polarization  switching  in  ferroelectric  ceramics  is  investigated. 
Under  bipolar  ac  electric  field  driving,  whereby  ferroelectric  ceramics  undergo  a  transition  between 
the  macroscopically  poled  and  depoled  states,  the  acoustic  emission  signal  is  large  in  the  poled  state 
and  small  in  the  depoled  state.  Under  unipolar  electric  field  driving,  whereby  a  prepoled  ferroelectric 
ceramic  sample  does  not  experience  a  macroscopic  polarization  reversal,  the  acoustic  emission  signal 
obtains  the  maximum  amplitude  at  an  intermediate  polarization  level.  Furthermore,  the  aeoustic  emis¬ 
sion  signal  decreases  as  the  induced  ferroelectric  polarization  decreases  after  ferroelectric  fatigue.  These 
phenomena  have  been  interpreted  by  attributing  the  AE  signal  to  the  elastic  energy  release  associated 
with  the  domain  switching  and  the  transmission  of  the  elastic  waves  via  a  electromechanical  coupling. 


INTRODUCTION 

Recently,  ferroelectric  switching  has  been  explored  for  a  number  of  potential  ap¬ 
plications.  Typical  examples  are  the  ferroelectric  thin  films  for  non-volatile  memory 
device.*  translucent  ferroelectric  ceramics  and  thin  films  for  optical  information 
storage,-  switchable  antiferroelectric  ceramics  for  high  strain  displacement  trans¬ 
ducer,’  and  soft  ferroelectric  ceramics  for  non-volatile  displacement  actuator.'’  Fer¬ 
roelectric  domain  switching  can  induce  a  longitudinal  strain  of  up  to  0.5%.  which 
can  cause  mechanical  failure  upon  electric  field  excitation.’  Consequently,  in-situ 
measurement  of  acoustic  emission  during  the  device  operation  has  been  proposed 
for  failure  prediction.’  In  an  early  work,  an  in-situ  measurement  has  indicated  that 
under  the  bipolar  field  driving  the  AE  signal  for  the  square  loop  PLZT(3/52/4‘^) 
was  most  pronounced  at  the  coercive  field  and  minimal  beyond  that.  The  authors’ 
attribute  this  to  the  pronounced  polarization  switch  around  the  coercive  field,  where 
the  macroscopic  depoling  took  place.  However,  for  PLZT(3/52/48)  ceramics  driven 
by  a  unipolar  ac  electric  field  and  for  relaxor  ferroelectric  PLZT(9/65/35)  driven 
by  a  bipolar  ac  field,  no  AE  signal  was  detected.  The  authors  attributed  the  missing 
of  AE  signal  to  the  absence  of  polarization  switching  in  these  two  cases. 

Acoustic  emission  signals  come  from  the  stress  waves  generated  during  the  events 
of  elastic  energy  release  within  a  solid  body.®  Domain  switching  and  crack  extension 
are  the  two  basic  mechanisms  of  the  acoustic  events  in  ferroelectric  ceramics  under 
applied  electrical  and/or  mechanical  fields.  The  energy  stored  in  the  stress  wave 
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must  be  transmitted  to  the  sensor/sample  interface  to  be  detected.  In  ferroelectric 
materials,  this  process  involves  electromechanical  coupling.  Consequently,  the  AE 
signals  detected  from  the  ferroelectric  ceramic  sample  should  be  a  function  of  the 
degree  of  poling.  This  behavior  has  not  been  revealed  in  a  previous  observation.^ 
Namely,  the  AE  signal  amplitude  is  maximal  at  the  polarization  switching  and 
minimal  when  the  poling  is  established.  It  is  the  intent  of  this  study  to  examine 
the  effect  of  electromechanical  coupling  on  the  detected  AE  signals  from  ferro¬ 
electric  ceramics  under  an  ac  electric  field  driving. 

PHYSICAL  BACKGROUND 

Acoustic  emissions  are  stress  waves  of  1  kHz  to  3  MHz  in  frequency,  that  are 
generated  inside  the  body  of  a  material  being  deformed  or  ruptured,  as  a  result  of 
elastic  energy  release  associated  with  the  processes.  The  stress  waves  are  attenuated 
within  the  body  usually  by  reflection  and  scattering.  By  attaching  a  PZT  sensor 
directly  to  the  specimen  at  deformation,  the  signals  can  be  correlated  to  the  energy 
release  associated  with  the  events.  The  initial  voltage  output  from  the  transducer, 
V,),  is  shown  to  be  proportional  to  the  square  root  of  the  energy  release, 
during  a  given  deformation  process,^ 

Vo  =  (1) 

where  >)/  is  a  proportionality  factor,  reflecting  both  the  preamplifier  characteristics 
and  the  attenuation  in  the  body  as  well  as  the  efficiency  of  coupling  between  the 
sensor  and  the  specimen.  For  materials  that  are  piezoelectric,  <1;  involves  an  ad¬ 
ditional  loss  due  to  the  electromechanical  coupling.  Specifically,  the  stress  waves 
would  be  converted,  through  direct  piezoelectric  effect,  to  electric  fields,  which 
would  then  be  converted  back  to  stress  waves  through  the  converse  piezoelectric 
effect.  Energy  transmission  is  involved  during  this  process,  as  governed  by  the 
electromechanical  coupling  factor  of  the  material.  It  is  known  that  ferroelectric 
ceramics  are  also  strong  piezoelectric  when  poled,  and  the  piezoelectricity  disap¬ 
pears  when  completely  depoled.  The  relaxor  ferroelectrics  are  piezoelectric  under 
a  DC  bias,  so  are  the  switchable  antiferroelectrics  when  the  applied  electric  field 
is  high  enough  to  cause  antiferroelectric-ferro-'iectric  phase  transition.  But  they 
are  not  piezoelectric  under  a  zero  bias.  Therefore,  the  electromechanical  atten¬ 
uation  of  the  AE  signals  within  a  ferroelectric  materials  is  a  function  of  the  po¬ 
larization  state. 

Under  an  ac  electric  field  driving,  the  ferroelectric  ceramics  experience  constantly 
poling  and  depoling  in  association  with  domain  switching.  A  large  mismatch  strain 
is  involved  in  the  process  due  to  the  high  anisotropic  mechanical  and  dielectric 
properties  along  different  crystallographic  directions,  leading  to  an  elastic  energy 
release  during  the  switching.  On  the  other  hand,  the  relaxor  ferroelectric  and 
antiferroelectric  do  not  exhibit  remanant  polarization  and  the  associated  domain 
structure  at  zero  bias.  The  relaxor  ferroelectrics  experience  a  microdomain  to 
macromain  phase  change  when  a  high  bias  is  applied.  Similarly,  the  antiferroelectric 
material  undergoes  an  antiferroelectric  to  ferroelectric  phase  change  under  a  high 
bias.  However  these  phase  change  processes  are  more  homogeneous  throughout 
the  body  and  involves  less  constraint  deformation.  As  a  result,  the  energy  associated 
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with  these  processes  is  less  than  that  involving  domain  switching  in  normal  ferro- 
electrics.  Another  important  acoustic  emission  event  involves  crack  propagation, 
both  from  the  macrocrack  and  microcracks.  When  microcrack  propagates  in  a 
ceramic  body,  the  energy  release  associated  with  cleavage  of  a  single  grain  of 
diameter  d  is,* 


where  £  is  the  Youngs  modulus  of  the  material  and  a  is  the  internal  stress  at  which 
fracture  takes  place. 

In  summary,  according  to  Equation  (1),  the  measured  signal  reflecting  Vq  will 
obtain  a  larger  value  for  an  event  involving  large  A£g  and  when  the  ferroelectric 
is  at  a  high  polarization  state. 


EXPERIMENTAL  PROCEDURES 

Ceramic  samples  studied  in  this  investigation  include  three  ferroelectric  ceramics: 
normal  ferroelectric,  switchable  antiferroelectric,  and  relaxor  ferroelectric  ce¬ 
ramics.  The  preparations  of  the  antiferroelectric  [Pbo  97Lao  o2(Zro.53Tio  ,2800  35)03], 
ferroelectric  [Pbo  ^sLao  osCZro  sgsTio  415)0  9903].  and  relaxor  ferroelectric  [0.9  Pb(Mgi,3 
Nb3y3)03:  0.1  PbTi03j  ceramic  samples  were  described  in  References  3,  4  and  8 
respectively.  Thin  specimens  were  cut  irom  each  of  the  three  sintered  Ceramic  body 
and  fine-ground  to  square  plates  of  dimensions  1  cm  x  1  cm  x  0.25  mm  and  1 
cm  X  0.5  cm  x  0.3  mm  for  the  antiferroelectric  and  the  ferroelectric  respectively, 
and  to  disc  plate  of  dimension  1.25  cm  x  1  mm  for  the  relaxor  ferroelectric 
ceramic.  The  thicknesses  of  the  samples  were  chosen  in  such  a  way  that  high  enough 
electric  field  can  be  applied  to  the  specimen  to  complete  ferroelectric  switching  for 
each  material.  The  electrode  was  put  on  the  specimen  either  by  a  metallic  silver 
sintering  technique  in  the  case  of  antiferroelectric  ceramic  or  by  a  gold  sputterig 
deposition  technique  for  normal  and  relaxor  ferroelectric  ceramics. 

The  dielectric  hysteresis  loops  of  these  ceramic  samples  were  obtained  using  a 
simple  Sawyer-Tower  circuit.  The  AE  signals  were  detected  using  an  acoustic 
emission  sensor  and  associated  amplification  system  (Acoustic  Emission  Inc.).  The 
samples  were  placed  on  the  sensor  with  one  major  surface  in  contact  with  the  sensor 
(Figure  1).  A  small  piece  of  alumina  was  placed  on  top  of  the  sample  to  ensure 
good  contact  between  the  sample  and  sensor.  In  addition,  a  thin  layer  of  organic 
coupling  agent  was  applied  in  between.  The  acoustic  emission  signal  and  the  P-E 
signals  from  the  Sawyer-Tower  system  were  simultaneously  acquired  using  an  IBM 
computer  with  a  16-bit  data  acquisition  board. 

The  PLZT  specimens  were  continuously  driven  by  a  1  Hz  electric  field  of  tri¬ 
angular  waveform  until  the  polarization  was  reduced  to  less  than  two  thirds  of  what 
was  obtained  from  the  fresh  sample,  a  process  referred  to  as  ferroelectric  fatigue.* 
Acoustic  emission  experiment  was  then  conducted  on  the  fatigue  specimens  without 
modifying  the  driving  condition. 


Polarization  (jjC/cm‘) 
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FIGURE  1  A  schematic  for  the  simultaneous  measurement  of  acoustic  emission  and  of  ferroelectric 
polarization. 


FIGURE  2  Dielectric  hysteresis  loops  (0.1  Hz)  of  three  typical  ferroelectric  types  of  ceramics:  a) 
Relaxor  ferroelectrics;  b)  Normal  ferroelectrics;  c)  Antiferroelectrics. 

RESULTS  AND  DISCUSSIONS 


The  P-E  hysteresis  loops  for  the  three  different  ferroelectric  ceramic  samples  are 
shown  in  Figure  2.  The  relaxor  ferroelectric,  normal  ferroelectric,  and  antiferro- 
electric  ceramic  samples  display  the  slim,  square  and  double  hysteresis  loops  re¬ 
spectively.  The  ferroelectric  polarization  and  AE  signal  in-situ  measurement  results 
for  the  three  different  ceramics  are  shown  in  Figure  3.  The  three  plots  have  one 
common  feature,  that  is.  the  AE  signal  amplitude  is  minimal  when  polarization  is 
being  reversed,  and  maximal  when  the  poling  is  established.  Clearly,  the  electro¬ 
mechanical  coupling  plays  an  important  role  in  governing  the  amplitude  of  the 
detected  AE  signal.  At  the  maximum  electromechanical  coupling  state  (poled 
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FIGURE  3  Acoustic  emission  signal  and  ferroelectric  polarization  versus  time  for  three  ferroelectric 
ceramics;  a)  Relaxor  ferroelectrics;  b)  Normal  ferroelectrics;  c)  Antiferroelectrics. 


State),  the  AE  signal  is  strong,  while  at  the  minimum  electromechanical  coupling 
state  (depoled  state)  the  AE  signal  is  weak. 

Figure  4  shows  the  result  of  the  in-situ  measurement  of  unipolar  field  induced 
polarization  and  the  corresponding  AE  signal  in  the  poled  5/58.5/41 .5  PLZT  ceramic 


124  W  PAN  and  H.  CAO 


0  12  3 


Time  (sec.) 

FIGURE  4  Acoustic  emission  signal  and  ferroelectric  polarization  versus  time  for  PLZT  (5/58.5/41 ,5) 
ceramic  under  a  unipolar  triangular  electric  field  driving. 


sample  (the  poling  field  is  in  the  same  polarity  as  the  unipolar  field).  Again  sub¬ 
stantial  AE  signals  were  observed  even  though  no  remnant  polarization  depoling 
was  experienced  by  the  sample  under  the  unipolar  field  driving.  In  this  case,  the 
acoustic  emission  signal  amplitude  is  more  uniformly  distributed  compared  with 
that  of  the  same  ceramic  under  the  bipolar  ac  field  driving.  This  is  because  the 
remnant  polarization,  which  is  the  origin  of  the  electromechanical  coupling,  is  not 
depoled  under  the  unipolar  ac  field  driving.  From  the  graph,  it  is  noted  that  the 
maximum  AE  signals  occurs  neither  at  the  maximum  induced  polarization  level 
nor  at  the  zero  polarization,  but  at  an  intermediate  polarization  level.  The  AE 
signal  variation  can  be  attributed  to  the  variation  of  the  domain  switching  under 
the  applied  electric  field.  Zhang  et  al.'*  demonstrated  that  in  poled  PZT  ceramics, 
the  irreversible  domain  wall  motion  can  be  excited  by  a  field  above  a  threshold 
field  which  is  smaller  than  the  macroscopic  coercive  field.  When  the  field  induced 
polarization  increases  from  zero  to  the  maximum,  the  polarizations  are  switched 
to  the  direction  of  the  applied  electric  field.  Near  the  zero  polarization  (zero  electric 
field),  the  applied  electric  field  is  too  small  to  switch  a  large  number  of  polarization 
vectors,  therefore  the  AE  signal  is  small.  Near  the  maximum  induced  polarization 
level,  on  the  other  hand,  most  of  the  switchable  polarization  vectors  are  already 
switched  to  the  direction  of  the  applied  electric  field,  the  polarizability  near  the 
maximum  polarization  level  come  largely  from  dielectric  permittivity  of  the  poled 
ceramics.  So  the  AE  signal  is  also  small.  The  polarization  switching  contributes 
more  to  the  polarizability  in  an  intermediate  {X>larization  level,  and  the  AE  signal 
is  thus  maximal  at  an  intermediate  polarization  level.  When  the  field  induced 
polarization  decreases  from  the  maximum  polarization  level  to  zero,  the  oriented 
polarization  vectors  drop  back.  That  occurs  before  the  applied  electric  field  de¬ 
creases  to  zero  because  these  polarization  vectors  are  not  a  part  of  the  remnant 
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FIGURE  5  The  acoustic  emission  signal  of  PL2T  (5/58-5/41 .5)  and  associated  polarization  versus  time 
for  a  partially  fatigued  fenoelecthc  (PLZT  5/58.5/41.5)  under  a  bipolar  triangular  wave  driving. 


polarization.  Therefore  the  maximum  AE  signal  also  occurs  at  an  intermediate 
polarization  level  as  the  driving  field  diminishes. 

The  effect  of  ferroelectric  fatigue  on  the  AE  signal  is  displayed  in  Figure  5  for 
PLZT  (5/58.5/41.5)  specimen.  When  the  ferroelectric  polarization  is  reduced  by 
fatigue,  the  AE  signal  decreases  accordingly.  This  is  consistent  with  above  obser¬ 
vations  on  three  ferroelectric  types  of  materials,  because  the  orientable  polariza¬ 
tion,  which  is  responsible  for  the  electromechanical  coupling  under  poled  state,  is 
reduced  after  fatigue. 

It  may  also  be  seen  from  Figure  3  that  the  AE  signal  for  the  antiferroelectric 
ceramic  sample  is  substantially  smaller  than  that  of  the  normal  ferroelectric  sample 
although  the  contact  area  of  the  antiferroelectric  ceramic  sample  is  larger  than  that 
of  the  normal  ferroelectric  ceramic  sample.  This  may  be  attributed  to  the  phase 
switching  in  antiferroelectric  ceramic  sample.  The  deformation  of  antiferroelectric 
ceramic  sample  under  the  applied  electric  field  is  more  homogeneous  inside  the 
body  and  consequently  fewer  AE  signals  are  generated. 


CONCLUDING  REMARKS 

Experiments  have  been  conducted  on  three  typical  ferroelectric  ceramics  to  in¬ 
vestigate  the  effects  of  electromechanical  coupling  on  the  acoustic  emission  under 
ac  electric  field  driving.  It  has  been  demonstrated  that  the  AE  signals  are  due  to 
the  events  of  domain  switching  and  the  signal  amplitude  is  strongly  dependent  upon 
the  polarization  state  in  the  body  as  a  result  of  the  electromechanical  coupling. 
Ferroelectrics  involving  domain  switching  generate  significantly  higher  AE  signal 
than  phase  change  ferroelectrics  (such  as  relaxor  ferroelectrics  and  antiferroelec- 
trics).  Under  the  same  macropolar  state,  the  AE  signal  amplitude  depends  on  the 
homogeneity  of  the  internal  strains. 
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ABSTRACT 


In  advanced  actuator  applications,  ferroelectric  ceramics  are  subjected  to  large 
alternating  electric  fields,  which  may  induce  cracking  and  mechanical  degradation 
phenomena.  This  airtide  describes  an  experimental  procedure  for  characterizing  crack 
extension  from  preexisting  flaws  in  ferroelectric  ceramics  subject  to  high  amplitude  but 
low  frequency  electric  fields.  Fracture  mechanics  concepts  are  applied  to  interpret  the 
observed  cracking.  Implications  for  the  mechanical  reliability  of  multilayered 
ferroelectric  actuators  are  discussed. 
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ABSTRACT 


For  ferroelectric  ceramics  used  either  as  sensors  or  as  actuators,  the  mechanical 
behavior  of  the  material  dictates  the  overall  electromechanical  response.  Furthermore, 
in  advanced  actuator  systems,  the  material  may  be  subject  to  high  stresses  and  large 
deformations,  leading  to  non-linearity  and  degradation  of  the  mechanical  performance. 
A  preliminary  experimental  investigation  has  been  undertaken  to  study  this  response. 
For  this  purpose,  polarization  and  strain  are  measured  continuously  on  a  compression 
cuboid.  Severe  non-linear  and  hysteretic  behavior  is  observed  in  many  cases.  An 
attempt  has  been  made  to  establish  the  constitutive  laws  applicable  to  this  non-linear 
deformation. 


KJS-Evans-27‘D«-OI*TA-Dfrm  Frict  91/07/10*4;36  PM*4/2/92 


2 


1.  INTRODUCTION 


Multilayer  ferroelectric  ceramics  are  used  as  actuators  and  sensors.  ^'2  Many 
advanced  applications  require  either  large  actuator  displacements  or  large  forces.  Such 
conditions  may  cause  non-linearity,  as  well  as  permanent  deformation,  of  both  the 
ferroelectric  and  the  metal  electrodes,  leading  to  degradation  of  the  performance  of  the 
actuator.3-5  Consequently,  it  is  essential  to  thoroughly  characterize  the  constitutive 
behavior  before  these  materials  can  be  successfully  integrated  into  systems.  Procedures 
for  establishing  the  constitutive  laws  are  developed  in  this  study.  The  non-linearities  are 
most  pronoimced  in  the  vicinity  of  stress  concentrations  that  occur  near  edges^  and  in 
association  with  flaws.^-S  The  latter  may  be  present  either  in  the  ferroelectric  or  at  the 
ferroelectric/electrode  interface.  Edges  as  well  as  flaws  in  the  ferroelectric  and  at 
interfaces  and  thus  afforded  particular  emphasis. 

Two  of  the  most  frequently  used  actuator/sensor  configurations  are  addressed: 
multilayer  and  annular  (Fig.  1).  The  multilayer  design  (Fig.  la)  has  been  used  for 
piezoelectrics  and  relaxor  ceramics^  as  well  as  capacitors.^  In  this  design,  the  stresses 
induced  by  the  actuator  displacements  are  aligned  with  the  polarization.  In  the  annular 
design  (Fig.  lb),  the  induced  stress  axis  is  normal  to  the  polarization  direction. 

Some  early  studies  indicated  that  poled  ferroelectrics  are  prone  to  non-linear 
deformation  upon  compression  along  the  poling  direction.^  As  a  result,  the  remnant 
polarization  is  partially  eliminated.  Conversely,  stresses  applied  perpendicular  to  the 
poling  direction  are  believed  to  stabilize  the  polarization.  The  non-linear  effects  are 
attributed  to  domain  switching,  induced  by  the  mechanical  loads.^0  A  more  detailed 
study  of  these  effects  is  presented  in  this  article.  An  attempt  is  also  made  to  rationalize 
the  observed  stress /strain  curves  with  respect  to  the  interplay  between  ferroelectric  and 
ferroelastic  domain  switching. 
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2.  EXPERIMENTAL  PROCEDURES 

2.1  Materials  and  Specimen  Preparation 

Ferroelectrics  considered  for  large  displacement  actuator  applications  are  of  three 
typesll-12  (Table  I):  i)  piezoelectrics;  ii)  ferroelectric  relaxors;  iii)  antiferroelectrics. 
Accordingly,  the  electromechanical  constitutive  laws  for  each  of  these  material  types  are 
explored.  Poled  piezoelectrics,  PZTs,  vere  obtained  from  a  commercial  sourcet:  PZT-A 
(Zr/Ti  =  52/48)  v\dth  hard  dopants  (Ba  and  Fe)  and  PZT-B  (Zr/Ti  =  53/47)  with  soft 
dopants  (Nb).  Relaxor  materials,  with  composition  0.9  PMN-0.1  PT,  were  prepared 
through  conventional  solid-state  sintering  and  by  hot-pressing.  The  antiferroelectric, 
with  a  composition  Pbo.97  La^  02  (Tio.64  Tio.09  Sr\)  27)  O3/  'vas  sintered  from  a  mixture  of 
reagent-grade  oxides. 

Short  compression  cuboids  with  nominal  dimensions  6x8x12  mm  were  cut  from 
bulk  ferroelectrics  and  all  non-electroded  faces  ground,  using  a  resin-bond  diamond 
wheel.  In  some  cases,  one  of  the  side  faces  was  polished  to  facilitate  observations  during 
testing.  For  the  relaxor  and  antiferroelectric  materials,  gold  electrodes  were  deposited 
onto  the  top  and  bottom  faces.  In  all  cases,  the  electrodes  were  connected  to  a  large 
capacitor  for  polarization  measurement. 

2.2  Testing  and  Measurements 

Uniaxial  compression  testing  of  brittle  materials  is  notoriously  difficult,  given  the 
requirement  that  the  load  be  applied  at  the  geometric  center  to  eliminate  bending. 
Nearly  perfect  alignment  is  essential  in  order  to  characterize  the  non-linear 
deformation.  Consequently,  a  special  loading  fixture  (Fig.  2)  has  been  designed  to 
facilitate  the  alignment  procedure  and  also  to  provide  a  means  for  recording  the  electric 
displacement.  Two  spherical  cones  were  incorporated  with  two  flat  alumina  blocks 

t  Channel  Industries,  Inc.,  Goleta,  California. 
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used  to  insulate  the  electrodes.  Specimens  were  carefully  emplaced  at  the  center 
between  the  alununa  blocks.  The  fixture  was  then  incorporated  into  a  servohydraulic 
machine  for  testing.  Final  adjustments  needed  to  minimize  bending  were  made  through 
in  situ  measurement  of  the  strain  differential  upon  exposing  the  specimen  to 
loading/unloading  cycles,  within  a  very  small  range  (<  5  MPa).  The  specimen  was 
connected  in  a  closed  loop  with  a  large  capacitor,  having  capacitance  five  orders  of 
magnitude  larger  than  the  specimen  (Fig.  2).  In  this  case,  the  voltage  change  across  the 
capacitor,  measured  by  a  voltage  follower,  provided  an  indication  of  the  charge  transfer 
from  the  specimen  to  the  capadtor.t  Consequently,  the  polarization  change,  AP,  could 
be  obtained  as 

AP  =  CqV  (1) 

where  Cq  is  the  capacitance  of  the  shunt  capacitor  and  V  is  the  voltage  measured  aaoss 
the  capacitor.  The  convention  adopted  here  is  that  depolarization  is  indicated  by  positive 
AP. 

Other  quantities  measured  simultaneously  with  the  voltage  include  the  load 
instantaneously  carried  by  the  specimen  obtained  from  the  load  cell  output  of  the 
testing  machine.  Also,  bi-directional  strain  gauges  were  emplaced  on  two  opposing  side 
faces  to  measure  the  deformation  along  both  the  loading  and  transverse  directions,  as 
well  as  monitor  possible  bending  caused  by  misalignment.  The  strain  component  along 
the  loading  direction  is  designated,  E^.  For  specimens  tested  with  the  loading 
perpendicular  to  the  polar  axis,  the  strain  along  that  axis  is  designated,  E^,  and  the  other 
normal  strain,  £y.  When  the  loading  is  parallel  to  polar  axis,  E,^  =  Ey.  Positive  quantities 
for  stress,  and  strain,  denote  compression,  unless  otherwise  indicated. 


f  It  is  emphasized  that  the  change  in  voltage  is  minimal. 
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3.  RESULTS 

3.1  Piezoelectric  PZT 

The  deformation  behavior  found  for  hard  and  soft  PZTs  subject  to  stress  applied 
parallel  to  the  poling  direction  has  similarities  and  important  differences  (Fig.  3).  Both 
materials  are  non-linear  when  the  stress  exceeds  ~  20  MPa,  and  exhibit  similar 
permanent  deformations  (-  0.5%)  after  exposure  to  stresses  above  500  MPa.  The  major 
difference  occurs  upon  unloading.  In  the  "soft"  material  (PZT-B),  the  non-linear 
deformation  is  nonrecoverable,  at  stress  levels  as  low  as  30  MPa  (Hg.  3a).  Conversely, 
the  "hard"  material,  P2rr-A,  recovers  most  of  the  non-linear  deformation,  with 
considerable  hysteresis  (Fig.  3b):  a  phenomenon  usually  referred  to  as  anelastic 
deformation.  However,  the  capacity  for  anelastic  strain  is  destroyed  at  stress  >  300  MPa. 

The  strain  measurements  indicate  that  the  volume  strains  (£z  +  Ex  +  6y)  are 
essentially  linear  for  both  material  (Fig.  4).  Consequently,  the  non-linear  deformation  is 
related  solely  to  the  deviatoric  components  of  stress,  a',  and  strain,  £'  (Fig.  5).  The 
instantaneous  Poisson's  ratio  (Fig.  6)  reflects  this  feature.  It  inaeases  initially,  reaches  a 
maximum,  decreases  upon  further  loading  and  finally,  approaches  a  value  characteristic 
of  incompressible  plastic  flow  (0.5).  The  deviatoric  stress /strain  curves  have  a  form 
(Fig.  5)  characteristic  of  that  found  for  phase  transformations,  consisting  of  i)  an  initial 
hardening  regime  wherein  the  hardening-rate  diminishes  as  the  stress  increases, 
followed  by  ii)  a  regime  of  accelerated  hardening,  as  a  saturation  plastic  strain  is 
approached.  Consequently,  the  stresses  a'  and  plastic  strain,  £'p,  can  be  described  by  a 
bifunctional  power  law  expression, 

=  (E;/e;)‘ +  (e^/e;)"  (2) 
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where  a  and  b  are  empirical  power  law  coefficients  (a  <  1  and  b  >  1),  is  a  reference 
stress  and  £  'q  a  reference  strain.  The  magnitudes  of  the  coefficients  for  the  hard  and  soft 
materials  are  summarized  in  Table  n. 

For  the  both  materials,  the  charge  on  the  capacitor  changes  during  loading, 
expressed  as  a  depolarization,  AP.  The  depolarization  curves  (Fig.  7)  for  each  material 
resemble  the  corresponding  stress/strain  curves  (Fig.  3),  suggesting  that  the  sources  of 
the  non-linear  deformation  are  also  responsible  for  the  polarization  changes.  The  total 
depolarization  measured  after  a  stress  exposure  of  ~  500  MPa  is  ~  80%  of  the  renmant 
polarization  induced  upon  p>oling. 

The  stress/strain  curves  obtained  when  the  stress  is  applied  normal  to  the  jxjling 
direction  are  qualitatively  the  same  (Fig.  8)  but  exhibit  more  limited  non-linearity.  For 
example,  the  permanent  strain  in  the  soft  material  (~  0.3%)  is  about  half  that  obtained 
upon  loading  along  the  poling  axis.  However,  the  transverse  strains  Ex  and  Ey  (Fig.  9) 
have  opposing  non-linearities,  as  a  result  of  the  anisotropy  caused  by  poling. 

The  depolarization  curves  exhibit  features  (Fig.  10)  dramatically  different  from  those 
evident  upon  loading  along  the  poling  axis.  Initially,  the  polarization  increases  as  the 
strain  increases,  but  when  the  strain  exceeds  ~  0.1%,  an  inversion  occurs  and  the 
polarization  begins  to  decrease  with  continued  straining. 

3.2  Rclaxor  PMN  PT 

The  uniaxial  compressive  stress-stredn  curves  obtained  for  PMN  PT  (Fig.  11)  indicate 
linear  elastic  behavior  up  to  an  axial  stress,  Gz  =•  500  MPa.  Thereafter,  some  non-linear 
deformation  appears.  Finally,  at  Gz  -  600  MPa,  fracture  occurs  and  the  specimen 
fragments.  Both  the  transverse  strain  and  volume  strain  (Fig.  12)  suggest  that  this  non¬ 
linearity  is  connected  with  microcracking.^^'^^  Notably,  the  volume  decreases  initially,  in 
accordance  with  linear  elasticity,  such  that  the  volume  strain  Ey  is  given  by. 
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e,  =  -(l-2v)a,/E 


(3) 


where  E  is  Young's  modulus  and  v  is  Poisson’s  ratio  (see  Table  I).  However,  as  non¬ 
linearity  begins,  dilatation  occurs  that  exceeds  the  elastic  shrinkage.  At  failure,  a  net  0.1% 
volume  expansion  is  evident.  Such  dilatations  are  typical  of  microcracking  effects  upon 
compression  loading  of  brittle  solids. 

There  were  no  polarization  changes  in  this  material,  consistent  with  the  PMN  PT 
being  mostly  paraelectric  (at  room  temperature)  and  having  no  remnant  polarization. 

3.3  Antiferroelectric  PLSnZT 

Appreciable  non-linear  deformation  is  observed  in  an  antiferroelectric  PLSnZT  at  a 
stress  slightly  higher  than  20  MPa*  (Fig.  13).  Some  of  this  deformation  is  permanent, 
manifest  in  the  unloading/ reloading  hysteresis.  Further  loading  (beyond 
Gz  =■  100  MPa)  causes  stiffening,  apparent  in  the  inaeased  tangent  modulus.  The 
hysteresis  loop  then  becomes  slimmer.  After  the  specimen  has  been  loaded  beyond 
~  500  MPa,  the  material  reverts  to  essentially  elastic  behavior.  The  total  permanent 
deformation  from  the  virgin  state  is  about  0.13%.  The  final  tangent  modulus  at 
saturation  is  -  160  GPa,  compared  with  an  initial  value  of  130  GPa.  The  compressive 
strength  is  ~  900  MPa  at  a  failure  strain  -  0.8%.  The  Poisson's  ratio  decreases  from  its 
initial  value  (V  -  0.35)  to  a  final  value  (~  0.25  at  strains  >  0.6%)  (Fig.  14),  typical  for 
oxide  ceramics. 

No  polarization  change  occurred,  consistent  with  the  overall  polarization  being  zero 
for  antiferroelectrics. 


*  No  orientation  distinction  is  necessary.  Antiferroelectric  materials  do  not  exhibit  remnant  polarization. 
They  are  isotropic  with  respect  to  mechanical  and  electrical  properties  at  zero  field  and  stress. 
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4.  NON-LINEAR  DEFORMATION 


Since  the  volume  strain  in  all  materials  app)ears  to  be  linear  and  non-hysteretic 
(unless  microcracking  occurs),  when  non-linearity  arises,  it  is  only  necessary  to 
understand  and  model  the  deviatoric  deformation:  as  in  metal  plasticity.  The 
incremental  deviatoric  strain  upon  loading  then  relates  to  the  corresponding  deviatoric 
stress,  through  a  constitutive  law  (Eqn.  2). 

The  mechanism  responsible  for  the  non-linear  deformation  in  the  piezoelectric  and 
antiferroelectric  materials  involves  90®  domain  switching,  having  the  features 
schematically  illustrated  in  Fig.  15.  The  as-processed  polycrystalline  material  has  a 
nominally  random  set  of  orientations  in  the  polarization  vector,  Pj.  An  applied  electric 
field  induces  both  180°  and  90®  domain  switching,  leading  to  a  preferred  alignment  of  Pj 
in  the  direction  of  the  field.  Residual  stress  also  develops  in  association  with  90®  domain 
switching.  Poling  of  FZTs  is  conducted  with  a  field,  E^.  Subsequent  stressing  may  be 
along  either  the  z  or  x.  For  domain  A  (Fig.  15),  oriented  along  z,  initially  Px  =  Py  =  0 
and  Pz  =  Pq-  This  type  of  domain  is  most  likely  to  switch  during  poling,t  above  a 
critical  field,  to  an  orientation  Py  =  Pz  =  0  and  Px  =  Pq-  The  new  orientation  is  stable, 
subject  to  subsequent  mechanical  stress  Gzz,  but  depoling  is  possible  when  a  axx  stress 
is  applied.  If  poling  did  not  cause  the  domain  to  switch  (because  the  field  is  too  small), 
switching  may  be  activated  by  a  subsequent  Gzz  stress,  whereupon  the  domain  would 
realign  along  either  Px  or  Py.  For  domain  C,  the  situation  is  similar  and  an  electric  field 
Ex  is  again  likely  to  induce  Py  ->  Px  switching.  Domain  B  is  most  illustrative  of  the 
interplay  between  the  electric  field  and  stress.  No  domain  switching  can  occur  subject  to 
Ex,  because  there  is  no  net  change  in  Px  up>on  switching  from  P^  to  P2:  moreover,  Px, 
decreases  if  switching  occurs  from  Pj  to  P3.  However,  upon  Gzz  stressing,  switching  from 


t  In  some  cases,  90**  domain  switching  upon  poling  may  be  inhibited  by  the  strain  energy  increase 
associated  with  the  residualization. 
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P|  to  P3  is  expected  when  the  normal  stress  (a|  =  7zz/2)  reaches  a  critical  value.  This 
event  produces  a  positive  strain  in  e  y  but  negative  strain  in  and  e^.  A  related  mode  of 
switching  may  occur  upon  Oxx  stressing. 

The  basic  insight  obtained  from  Fig.  15  allows  the  major  phenomena  found  by 
experiment  to  be  rationalized.  For  the  poled  PZT  materials,  tested  with  the  mechanical 
loading  axis  coincident  with  the  poling  direction,  the  mechanical  field  causes  switching 
in  a  sense  opposite  to  that  caused  by  poling  resulting  in  systematic  depoling.  Once  the 
depoling  is  complete,  further  switching  does  not  occvir.  However,  this  mechanically 
depoled  state  must  be  contrasted  and  distinguished  from  the  original  (unpoled)  state,  in 
that  the  spaciously  random  polarization  vectors  in  the  virgin  material  are  now  in  a 
plane  normal  to  the  depoling  stress.  "Complete  depoling"  results  in  the  saturation 
behavior  found  at  plastic  strains  >  0.5%  (Fig.  3).  Thereafter,  the  material  behaves  as  an 
elastic  body  until  microcracking  commences  prior  to  mechanical  failure. 

Loading  normal  to  the  poling  axis  (usually  thought  to  be  linear)  causes  a  non-linear 
response  best  illustrated  by  orientation  B  in  Fig.  15.  Some  of  the  90°  domains  are 
unfavorably  oriented  for  switching  upon  poling,  but  are  in  a  preferred  orientation  for 
switching  upon  loading.  The  important  consequence  is  the  relationship  found  between 
the  polarization  change,  AP,  and  the  strain  (Fig.  10).  The  initial  portion  of  these 
curves  is  the  normal,  linear  piezoelectric  effect.  However,  type  B  domain  switching  at 
larger  stresses  causes  depolarization  (-  APx)  and  an  expansion  along  the  x-axis,  leading 
to  the  observed  inversion  in  the  polarization  curve. 

The  parameters  that  characterize  the  non-linear  deviatoric  stress/strain  (Eqn.  2), 
must  depend  on  factors  such  as  the  domain  wall  mobility,  the  residual  stresses,  etc. 

More  detailed  analysis  is  needed  to  develop  a  general  yield  criterion  and  a  flow  law.  For 
this  purpose,  experiments  involving  multiaxial  stresses  are  necessary.  Furthermore,  the 
differences  found  between  the  "hard"  and  "soft"  PZTs  should  be  consistent  with  the 
yield  criterion  and  flow  law.  One  key  feature  concerns  the  residual  stresses,  which  are 
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expected  to  be  larger  in  the  hard  material,  causing  the  domains  to  switch  back  when  the 
applied  load  is  removed. 

For  the  antiferroelectric  material,  the  permanent  strain  is  associated  with  ferroelastic 
domain  switching,  with  characteristics  similar  to  those  described  above  for  PZT. 
However,  it  is  noted  that  electric  fields  cause  the  antiferroelectrics  to  ferroelectric  phase 
transition. 


5.  CONCLUSION 

Mechanical  non-linearity  found  in  piezoelectric  and  antiferroelectric  materials  has 
been  attributed  to  ferroelastic  domain  switching,  with  a  yield  criterion  :  jntrolled  by  a 
deviatoric  effective  stress.  The  flow  rule  has  not  been  established.  However,  curve 
fitting  suggests  a  power  law  at  small  strains,  similar  to  metal  plasticity.  At  larger  strain, 
accelerated  hardening  occurs  and  the  inelastic  strain  approaches  saturation.  In  addition, 
strong  interconnections  have  been  established  between  ferroelectric  and  ferroelastic 
deformation.  Notably,  electrically  "hard"  PZT  can  sustain  a  large  mechanical  load 
without  significant  permanent  deformation,  whereas  a  "soft"  PZT  is  susceptible  to 
deformation  at  small  stresses,  even  though  the  "yield"  stresses  are  similar  for  the  two 
materials.  Consequently,  for  actuator  applications,  hard  PZT  can  be  used  at  higher 
stress  levels,  without  depoling. 

By  contrast,  relaxor  ferroelectric  ceramics  are  capable  of  sustaining  large  mechanical 
stresses  without  either  mechanical  degradation  or  non-linearity.  Furthermore,  the 
strain-polarization  curves  have  negligible  hysteresis.  These  characteristics  render  the 
material  a  good  candidate  for  large  displacement,  large  force  actuators.  In  addition, 
among  the  materials  investigated,  this  material  has  the  highest  modulus,  making  it 
possible  to  produce  appreciable  actuation  displacements,  while  sustaining  relatively 
large  mechanical  stress. 
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TABLE  I 


Materials  and  Properties 


^Dlelactrlc 

Constant 

Coeciva 

Field 

Sictti  rA.te/ 

Polarization 

Elastic 

Constants 

Piezoelectric 
Constants  (KT'^mAO 

Material 

Pr 

(C/m^) 

1 

4 

(GPa) 

d33 

d3i 

di5 

(Hard) 

PZT-A<» 

1300 

2.f 

0.35 

65 

82 

25 

300 

-135 

525 

(Soft) 

PZT-B» 

1750 

/  ^ 

0.35 

52 

64 

20 

400 

-185 

625 

PLSnZT 

^00 

5.1 

0.31 

140 

PMN*PT 

20.000 

B 

255 

^  K  is  measured  at  room  temperature  and  10  kHz. 
^  Channel  Industries,  Inc.,  Goleta,  California. 
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TABLE  11 


Constants  in  Deviatoric  Constitutive  Law  (Eqn.  2) 


Exponent 

a 

Exponent 

b 

Reference 

Stress 

O'o  <MPa) 

Reference 

Strain 

e'o(%) 

(Hard)  PZT-A 

0.54 

9 

100 

0.62 

(Soft)  PZT-B 

0.7 

12 

75 

0.65 

PlSnZT 

0.65 

9 

100 

0.2 
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FIGURE  CAPTIONS 

Fig.  1.  Actuator  configuratioiis:  a)  multilayer,  b)  annular. 

Fig.  2.  Compression  loading  fixture  with  associated  measurement  circuits  and 
orientation  designation. 

Fig.  3.  Stress/ strain  curves  from  hard  and  soft  PZTs  subject  to  stress  along  polar  axis. 

Fig.  4.  Change  in  volume  strain  with  axial  stress  for  hard  and  soft  PZTs. 

Fig.  5.  The  deviator  plastic  strain,  as  it  relates  to  deviatoric  stress,  is  solely  responsible 
to  the  non-linear  deformation. 

Fig.  6.  Change  in  instantaneous  Poisson's  ratio  with  axial  strain  for  soft  PZT. 

Fig.  7.  Electrical  depolarization,  AP,  for  hard  and  soft  PZTs  as  a  function  of  the 
applied  stress. 

Fig.  8.  Stress/strain  curves  for  hard  and  soft  PZTs  when  the  stress  is  applied  normal 
to  the  polar  axis. 

Fig.  19.  Transverse  strain  characteristics  for  both  PZTs  as  a  function  of  applied  normal 
stress. 

Fig.  10.  Polarization  changes  in  the  polar  axis  caused  by  stress  applied  normal  to  the 
polar  axis. 

Fig.  11.  Uniaxial  stress/strain  curves  for  PMN  •PT  relaxor  material. 

Fig.  12.  Effects  of  stress  on  a)  transverse  and  b)  volume  strain  for  PMN»PT. 

Fig.  13.  Uniaxial  stress/strain  curve  for  the  antiferroelectric  PLSnCT. 

Fig.  14.  Change  in  the  Poisson's  ratio  as  a  function  of  stress  for  PLSnZT. 

Fig.  15.  Schematic  of  deformation  by  domain  switching. 
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Abstract 


In  advanced  actuator  systems,  ferroelectric  ceramics  are  usually  made  into 
multilayered  stacks  to  obtain  optimum  performance.  Piis  is  done  by  interlacing 
noble  metals  (such  as  Pt )  with  ceramic  green  sheet  and  cofiring.  The  present  study 
proposes  an  alternative,  more  economic  electrode.  Copper  has  been  successfully 
bonded  to  several  PZT-based  ferroelectric  ceramics  in  a  multilayer  configuration  by 
diffusion  bonding.  The  microstructure  and  the  relevant  physical  properties  of  the 
interfaces  have  been  investigated  using  microscopy  as  well  as  mechanical  and 
electrical  characterizations.  These  studies  indicate  that  no  reaction  product  forms  at 
the  interface  and,  furthermore,  the  interface  exhibits  high  resistance  to  decohesion, 
making  Cu  a  promising  alternative  electrode  material. 
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1.  INTRODUCTION 


The  application  of  ferroelectric  ceramics  as  actuators  and  sensors  is  reliant  on 
multilayer  concepts,  wherein  electrodes  alternate  with  the  ferroelectric. ^'2  while  the 
electrodes  are  often  regarded  as  being  secondary  to  the  actuator/ sensor  system,  the 
ultimate  performance  of  the  system  is  strongly  influenced  by  the  electrical  and 
mechanical  characteristics  of  the  electrode  itself  and,  usually,  the  electrode/ 
ferroelectric  interface.^  The  interfacial  zone  has  particular  ramifications  on 
performance,  through  a  multiplicity  of  potential  degradation  effects:  fatigue, 
decohesion  and  charge  injection.  There  has  been  minimal  basic  research  on 
interfaces  between  metals  and  complex  oxides  (other  than  binary  systems,  AI2O3, 
NiO,  etc.)^.  Consequently,  there  are  no  rules  on  such  issues  as  reaction  products, 
decohesion,  interface  resistance,  etc.  Furthermore,  there  have  been  few  systematic 
studies  of  ferroelectric /electrode  interfaces  and,  as  yet,  there  are  no  scientific 
governing  principles. 

Common  electrode  systems  for  PZT  sensors  and  BaTi03  capacitors  are  typically 
noble  metals  (such  as  Pt  and  Ag  ).  These  materials  are  expensive  and  may  be 
susceptible  to  degradation  effects^,  such  as  debonding.  Furthermore,  Ag  is  subject  to 
migration  effects  that  lead  to  prematvire  breakdown.  Consequently,  a  Pd/Pt  alloy, 
that  is  apparently  resistant  to  migration,  is  often  used  for  cofired  multilayer 
actuators.  The  present  article  represents  an  initial  attempt  at  studying  a  rational 
interface  between  a  ferroelectric  and  an  inexpensive  metal,  with  implications  for 
multilayer  systems. 

For  this  purpose,  bonds  consisting  of  alternate  layers  of  Cu  and  two  ferroelectric 
ceramics  have  been  produced  by  diffusion  bonding.  The  microstructure  of  the 
interface  has  been  characterized  using  a  combination  of  electron  microscopy 
techniques.  The  interfacial  mechanical  properties  have  also  been  assessed,  using  an 
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indentation  technique.  Furthermore,  the  Polarization /Field  (P/E)  behavior  has  been 
characterized  in  order  to  validate  that  Cu  can  be  used  as  an  economic  alternative 
electrode,  without  adverse  effects  on  the  ferroelectric  and  electromechanical 
response. 

2.  EXPERIMENTAL 

2.1  Materials 

The  oxide  ferroelectrics  used  for  this  study  include  a  hard  PZT^  and  an 
antiferroelectric  PLSnZT,  with  a  chemical  composition  Pb  97La  o2(Zr  66Ti  |iSn  23)03. 
Thin  plates  having  dimensions  5x8  mm  and  2  mm  thickness  were  cut  from  bulk 
sintered  material  and  ground  using  a  resin-bonded  diamond  wheel  to  mininuze 
machining  damage.  High  purity  (99.98%)  Cu  foil,  25  |J.m  in  thickness,  was 
interspersed  to  form  a  multilayer  (Fig.  1).  Bonding  was  conducted  in  a  vacuum  at 
1050‘’C  for  10  hrs.  A  small  pressiue  of  about  5  MPa  was  maintained  throughout. 
Under  these  conditions,  lead  loss  is  minimal.  A  slow  cooling  rate  of  about  3°C/min 
was  used  to  alleviate  residual  stress.  The  bonded  specimen  was  sectioned  normal  to 
the  interface  and  polished. 

2.2  Characterization 

Thin  sections  for  transmission  electron  microscopy  (TEM)  were  prepared  using 
the  off-center  dimpling  technique  discussed  elsewhere.^  Discs,  3  mm  in  diameter, 
with  the  interface  region  about  0.5  mm  from  the  center  of  the  dimple  were  ion- 
milled  at  4  kV,  15  ‘  beam  angle.  Semi-circular  quartz  shields  were  used  to  protect 
one  side  of  the  interface  from  the  ion  beams.  The  resulting  samples  showed  a  large 
hole  adjacent  to  the  interface  in  the  ceramic  (about  0.8  mm  diameter)  with  interface 


+  PZT  C5500,  Chanel  Industries  ,  Inc.  Goleta,  California 
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regions  sufficiently  thin  for  TEM  observations  at  400  kV.  Electron  diffraction 
measurements,  energy  dispersive  X-ray  analysis  (EDS),  high  resolution  imaging 
(HREM)  and  parallel  electron  energy  loss  spectroscopic  measurements  (PEELS)  were 
all  carried  out  on  a  JOEL  4000  FX  side  entry  microscope. 

Bonded  samples  with  Cu  electrodes  were  incorporated  into  a  Sawyer-Tower 
circuit  to  measure  the  polarization  (P/E)  behavior.  For  comparison,  the  P-E  curves 
of  the  ferroelectrics  were  also  measured  by  using  thin  slice  samples  with  sputtered 
gold  electrodes.  The  capacitance  of  the  multilayer  was  measured  by  an  LRC  meter. 
Dielectric  breakdown  measurements  were  also  performed,  by  applying  a  high  DC 
voltage  across  the  specimen. 

2.3  Indentation 

The  fracture  resistance  of  the  interface  between  the  metal  electrode  and  the 
ferroelectric  was  assessed  using  an  indentation  technique.  A  Vicker's  indent  was 
emplaced  into  the  ceramic,  adjacent  to  the  interface,  aligned  such  that  one  radial 
crack  had  normal  incidence  with  the  interface  (Fig.  2).  Indentation  loads  between  5 
and  30N  were  used.  The  consequent  deformation  of  the  metal  and  slip  along  the 
interface,  as  well  as  the  occurrence  of  interface  decohesion,  were  inspected  in  a 
scanning  electron  microscope,  after  the  specimen  had  been  coated  with  gold  by 
plasma  sputtering. 

3.  RESULTS 

3.1  Interface  Structure 

i)  Cu  /  PLSnZT 

A  low  magnification  TEM  view  of  Cu  /  PLSnZT  (Fig.  3a)  indicates  a  relatively 
planar  interface.  Faceted  regions  were  identified  at  several  locations  (Fig.  3b). 
Diffraction  and  EDS  established  that  these  regions  were  pure  Cu,  with  no  evident 
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interdiffusion  with  the  ferroelectric.  The  Cu  grain  size  was  substantially  larger  than 
the  ferroelectric  grain  size.  PEELS  measurennents  confirnned  the  absence  of  O  in  the 
Cu  layer,  although  in  some  Cu  grains  close  to  the  interface,  a  weak  N-edge  was 
visible.  Observations  at  high  magnification  were  unable  to  identify  a  discrete  Cu 
oxide  layer  at  the  interface,  but  revealed  the  presence  of  Zr02  precipitates  (Fig.  4a, 
arrowed  regions,  and  b).  The  shape  of  the  precipitate  suggests  a  cubic  crystal 
structure.  The  composition  of  the  precipitate  was  established  by  means  of  combined 
EDS  and  PEELS  analysis.  The  average  precipitate  size  was  about  50  nm  and  all 
particles  were  observed  to  protrude  into  the  metal  layer. 

Electron  diffraction  analysis  of  the  ferroelectric  perovskite  [001]  and  [110]  zone 
axis  orientations^/®  revealed  the  presence  of  a  long  period  superlattice  modulation 
of  an  (as  yet)  undetermined  origin.  A  complete  structure  determination  of  this 
phase  by  means  of  electron.  X-ray  and  neutron  diffraction  is  currently  underway. 
High  defect  concentrations  (dislocations  and  domain  walls)  were  observed  in  the 
ferroelectric  close  to  the  interface  (Fig.  4a). 

ii)  Cu  /  PZT 

The  interface  in  the  cross-section  Cu  /  PZT  samples  is  not  as  planar  as  that  in 
the  Cu  /  PLSnZT  samples  (Fig.  5a).  EDS  analysis  of  the  metal  again  revealed  no 
interdiffusion  of  cations  from  the  ferroelectric  into  the  Cu.  Several  precipitates  occur 
along  the  interface  having  sizes  in  the  range  40  to  100  nm.  Electron  diffraction 
analysis  combined  with  EDS  identified  these  particles  as  cubic  zirconia.  Again,  a 
high  defect  concentration  was  observed  in  the  ferroelectric  close  to  the  interface, 
consisting  of  both  the  dislocations  and  domain  cells  (Fig.  5b).  The  domain  walls  are 
identified  by  the  fringe  contrast  with  features  being  asymmetric  in  bright  field  and 
symmetric  in  dcirk  field  imaging  mode.  In  addition,  well-developed  domains  were 
present  in  all  ferroelectric  grains. 
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3.2  Measurements 


Capacitance  measurements,  performed  on  the  multilayer,  gave  values  for  the 
dielectric  constants,  K  =  1460.  This  value  is  identical  to  that  measured  on  the 
material  with  sputtered-gold  electrodes,  prior  to  diffusion  bonding.  It  also  coincides 
with  the  value  quoted  by  the  manufacturer  The  ferroelectric  hysteresis  loop,  and 
the  dielectric  strength  measured  on  the  bonded  trilayer  system,  were  also  identical  to 
values  obtained  before  bonding. 

Cracks  introduced  within  the  ferroelectric  by  indentation  that  extend  to  the 
electrode  interface  are  invariably  arrested  at  that  interface  (Fig.  6)  .There  is  no 
evidence  of  interface  debonding.  Furthermore,  the  plastic  zone  developed  in  Cu  is 
similar  to  those  developed  in  an  Al/ Alumina  systemlO;  a  well-characterized  strong 
interface  between  a  metal  and  an  oxide.  Typically,  the  cracks  blunt  at  the  interface, 
indicative  of  a  'strong'  bond  and  extensive  yielding  in  the  metal.  Such  strong 
interface  bonding  has  also  been  noted  between  Cu  and  certain  binary  oxides.^ 


4.  ANALYSIS 

Creation  of  a  thin  foil  section  in  preparing  the  TEM  samples  introduces 
surfaces  that  act  as  edge  sites,  which  concentrate  the  thermal  expansion  mismatch 
stress!  1.  Yielding  occurs  from  this  edge  if  the  equivalent  stress,  Oe,  exceeds  the  yield 
strength  of  either  constituent  material.  Dimensional  arguments  suggest  that  the  foil 
thickness,  t,  relative  to  the  thickness  of  the  Cu  layer,  h,  controls  plastic  zone 
evolution,  as  thinning  proceeds  (Fig.  7a).  Consequently,  two  typical  t/h  ratios  are 
analyzed:  1  and  0.1.  The  latter  represents  a  limiting  case,  wherein  the  bond  can  be 
regarded  as  a  bi-material  interface  separating  semi-infinite  Cu  and  ferroelectric 
materials. 
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A  finite  element  procedure  (ABAQUS)  has  been  used  to  analyze  the  plastic 
strain  contours  at  the  interface.  The  elastic-plastic  behavior,  assumed  for  Cu,  is 
governed  by  a  von  Mises  yield  condition,  with  power  law  hardening  (Table  I).  For 
the  ferroelectric  material,  a  two-term  power  law  is  used  to  model  the  ferroelectric 
deformation,  associated  with  domain  switching  under  stress  ^2.  Constitutive 
parameters  used  for  this  calculation  are  summarized  in  Table  I.  Geometric 
symmetry  allows  the  calculations  to  be  performed  using  a  quarter  of  the  cross- 
section(Fig.  7).  The  plastic  strain  contours  are  computed  (for  both  t/h 
configurations),  as  the  specimen  is  cooled  from  the  processing  temperature.  It  is 
assumed  that  creep  relaxation  occurs  at  high  temperatures,  such  that  residual 
stresses  do  not  begin  to  develop  above  T  «  0.7  T^^,  with  T^  being  the  melting  point  of 
the  copper.  The  results  (Fig.  8)  indicate  a  relatively  small  effect  of  t/h.  In  both  cases, 
there  is  a  much  larger  plastic  zone  in  the  ferroelectric  than  in  the  Cu.  The  plastic 
strain  in  the  PZT  near  the  interface  approaches  0.6%  (  the  permanent  saturation 
strain  exhibited  by  PZT  ^2).  The  features  exhibited  by  the  plastic  zone  are  consistent 
witn  the  high  dislocation  density  and  the  domain  wall  configurations  observed  by 
TEM  (Fig.  4,  5). 


5.  Concluding  Remarks 

Two  of  the  present  observations  indicate  that  Cu  has  desirable  characteristics  as 
an  electrode  for  a  ferroelectric,  multilayer  actuator,  i)  There  are  no  reaction  products 
at  temperatures  needed  for  diffusion  bonding,  (except  for  occasional,  small 
precipitates  of  Zr02).  h)  The  interface  is  strongly  bonded  and  resistant  to  decohesion. 

A  possible  detriment  is  the  occurrence  of  plastic  deformation  in  the 
ferroelectric,  near  the  interface.  Such  deformation  is  «m  edge  effect,  associated  with 
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the  thermal  expansion  mismatch,  and  arises  because  of  the  low  'yield  strength'  of 
the  ferroelectric.  The  concern  is  that  superposed  effects  of  electric  field,  during 
operation,  may  result  in  fatigue  damage  that  spreads  from  the  edge^3.  Further  study 
is  deemed  necessary  to  elucidate  this  effect. 
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Table  I  Contitutive  Properties  of  the  Cu  /  PZT  System 


Material 

Cu 

PZT 

Young's  modulus,  E(GPa) 

130 

80 

Poisson's  ratio,  v 

0.343 

0.4 

Thermal  expansion 

1.8  (T>Tc) 

coefficient,  a  (lO'^K'l) 

16.6 

9.0  (T<Tc) 

Yield  strength,  Oq  (MPa) 

60 

a 

Hardening  exponent,  n 

0.32 

Reference  stress,  o'q  (MPa)§ 

75 

Reference  strain,  e'^  ^ 

0.65 

Exponent  a  ^ 

0.7 

Exponent  b  ^ 

12 

The  power  hardening  law  for  Cu  relates  the  true  stress,  a,  with  the  true  plastic  strain,  e; 
a  =  Oo(e/eo)" 

The  two-term  power  law  for  PZT  satisfies  the  following  expression 

o7cJ'o=(e7e'or  +  (e7e'o)‘’ 


where  o'  and  c'  are  the  deviatoric  stress  and  plastic  strain  ,  respectively. 


KJS-Evans-5«S*TA-Str  Prp  91/07/09«13:08*7/19/92 


10 


REFERENCES 


[1]  S.R.  Winzer,  N.  Shanker  and  A.P.  Ritter,  "Designing  Cofired  Multilayered 
Electrostrictive  Actuators  For  Reliability,"  /.  Am.  Ceram.  Soc.,  72,  2246-57 
(1990). 

[2]  L.  E.  Cross,  "Piezoelectric  and  Electrostrictive  Sensors  and  Actuators  for 
Adaptive  Structures  and  Smart  Materials."  ONR  Annual  Repart,  Penn.  State 
University,  (1990). 

[3]  W.J.  Menz  and  R.  Anderson,  AT&T  Bell  Laboratory  Report.  (1978). 

[4]  I.E.  Reimanis,  B.  J.  Dalgleish  and  A.G.  Evans,  "The  Fracture  Resistance  of  A 
Model  Metal/Ceramic  Interface,"  Acta  metall.  39,  3133-3141, 1991. 

[5]  M.  Ruble,  A.  G.  Evans,  M.  F.  Ashby  and  J.  P.  Hirth,  Metal  Ceramic  Interfaces, 
Acta,  Scripta  Metallurgica  Proc.  Series,  Vol.  4  (1990). 

[6]  M.  DeGraef,  B.J.  Dalgleish,  M.R.  Turner  and  A.G.  Evans,  "Interfaces  Between 
Blumina  and  Platinum:  Structure  and  Fracture  Resistance,"  Acta  Metall. 
Mater,  in  press,  (1992). 

[7]  Y.-J.  Chang,  Y.-J.  Lian  and  Y.-L.  Wang,  "One  dimensional  regular  arrays  of 
antiphase  domain  boundaries  in  anti-ferroelectric  Sn-substituted  lead 
zirconate  titanate  (PZT)  ceramics,"  Appl.  Phys.,  A36,  221  (1985). 

[8]  M.  De  Graef,  J.  Speck,  D.R.  Clarke  and  D.  Dimos,  "Electron  microscopy  study 
of  domain  boundaries  in  releixor  ferroelectrics,"  submitted  to  Symposium  I, 
Ferroelectric  Thin  Films,  Materials  Research  Society  Fall  Meeting  (Boston 
1991). 


KJS-Evans-5*S*TA-Str  Prp  9 1/07/09*  13:08* 7/1 9/92 


11 


[9]  Catalogue,  Channel  Industries,  Inc. 


[10]  H.  C.  Cao,  M.  D.  Thouless  and  A.  G.  Evans,  "Residual  Stresses  and  Cracking  in 
Brittle  Solids  Boded  with  a  Thin  Ductile  Layer."  Acta  metall.  36  [8],  2037-2046, 
(1988). 

[11]  H.  C.  Cao,  "A  Model  for  Crack  Extension  in  Metal  Reinforced  Brittle  Solids," 
Advanced  Composites  Materials,  Ceramic  Transactions  (Ed.  M.  Sacks)  Vol.  42, 
pp.  485-492,  (1991). 

[12]  H.  C.  Cao  and  A.  G.  Evans,  "Nonlinear  (Deformation  of  Ferroelectric 
Ceramics."  submitted  to  /.  Amer.  Ceram.  Soc.,  (1992). 

[13]  Q.  Y.  Jiang,  L.  E.  Cross,  et  al,  to  be  published. 

[14]  J.  H.  Hollomon,  "Tensile  Deformation,"  Trans.  AIME,  162,  p.  268,  (1945). 


KJS-Evans-5*S*TA-Str  Prp  91/07/09*  13:08*7/1 9/92 


12 


FIGURE  CAPTIONS 


Fig.  1.  Schematic  showing  the  layout  of  Cu/PZT  multilayered  stack  for  diffusion 
bonding 

Fig.  2.  A  schematic  illustrates  the  indentation  in  a  ferroelectric  ceramic  close  to  an 
interface  and  the  associated  plastic  zone  in  the  electrode. 

Fig.  3.  a)  Low  magnification  TEM  micrograph  indicates  a  planar  interface 
between  Cu  and  PZT  (arrows  indicate  locations  of  the  interface), 
b)  Faceted  Cu  regions  near  the  Cu/PLSnZT  interface.  Some  precipitates 
can  also  be  observed  near  the  interface. 

Fig.  4.  a)  50  nm  Zr02  precipitates  (Arrowed)  at  the  Cu/PLSnZT  interface. 

b)  A  cube-shaped  Zr02  precipitate  at  the  interface,  protruding  into  the 
metal  layer. 

Fig.  5.  a)  Interface  (arrowed)  between  Cu  an  d  PZT.  Ferroelectric  domain  walls 
are  present  in  all  grains. 

b)  High  defect  concentration  in  the  PZT  close  to  the  interface. 

Fig.  6.  Indentation  crack  arrests  at  the  electrode/ ferroelectric  interface  and 
blunting  takes  place.  No  interface  decohesion  is  observed  (interfaces 
marked  by  arrows). 

Fig.  7.  Geometry  of  thin  foils  used  for  FEM  analysis  and  a  typical  finite  element 
mesh  for  t/h=l. 
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Fig.  8.  Plastic  strain  contours  in  a  Cu/PZT  bond,  caused  by  thermal  mismatch 
straining:  a)  t/h  =  1;  b)  t/h  =  0.1. 
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Non-Linear  Constitutive  Properties  of 
Piezoelectric  Ceramics 

H.C.  Cao  and  A.G.  Evans 

Materials  Department  University  of  California,  Santa  Barbara,  California  93106 


ABSTRACT 

Experiments  have  been  conducted  to  study  the  nonlinear  constitutive 
properties  of  ferroelectric  ceramics  for  large  strain  operation  applications,  where  the 
mechanical  behavior  of  the  material  dictates  the  overall  electromechanical  response 
of  the  devices.  For  this  purpese,  pelarization  and  strain  are  measured  continuously 
on  a  compression  cuboid,  ^vere  non-linear  and  hysteretic  behavior  is  observed  in 
these  materials.  An  empirical  formulation  has  been  proposed  to  model  the 
constitutive  laws  applicable  to  this  non-linear  deformation.  Physical  mechanism  is 
provided  ,  that  is  attributed  to  this  constitutive  law. 

1.  INTRODUCTION 

Multilayer  ferroelectric  ceramics  are  used  as  actuators  and  sensors.’-^  Many 
advanced  applications  require  either  large  actuator  displacements  or  large  forces. 
Such  conditions  may  cause  non-linearity,  as  well  as  p^ermanent  deformation,  of  both 
the  ferroelectric  and  the  metal  electrodes,  leading  to  degradation  of  the  performance 
of  the  actuator.3‘5  The  non-linearities  are  most  pronounced  in  the  vicinity  of  stress 
concentrations  that  occur  near  edges^  and  in  association  with  flaws.^'®  The  latter 
may  be  present  either  in  the  ferroelectric  or  at  the  ferroelectric /electrode  interface. 
Consequently,  it  is  essential  to  thoroughly  characterize  the  constitutive  behavior, 
before  these  materials  can  be  successfully  integrated  into  system.*;.  Some  of  the 
procedures  necessary  for  establishing  the  constitutive  laws  are  develojjed  in  a  series 
of  studies. 

Two  of  the  most  frequently  used  actuator/sensor  configxirations  are  addressed: 
multilayer  and  annular,  ^th  designs  have  been  used  for  piezoelectric  actuators.^  In 
the  multilayer  design,  the  stresses  induced  by  the  actuator  displacements  are  aligned 
with  the  polarization.  In  the  annular  design,  the  induced  stress  axis  is  normal  to  the 
polarization  direction. 

Some  early  studies  indicated  that  poled  ferr.  ■electrics  are  prone  to  non-linear 
deformation  upon  compression  along  the  poling  direction.^  As  a  result,  the  remnant 
polarization  is  partially  eliminated.  Conversely,  stresses  applied  perpendicular  to  the 
poling  direction  are  believed  to  stabilize  the  polarization.  TTie  non-linear  effects  are 
attributed  to  domain  switching,  induced  by  the  mechanical  loads.^0  a  more  detailed 
study  of  these  effects  is  presented  in  this  article.  An  attempt  is  also  made  to 
rationalize  the  observed  stress/strain  curves  with  respoct  to  the  interplay  between 
ferroelectric  and  ferroelastic  domain  switching. 

2.  EXPERIMENTAL  PROCEDURES 

2.1  Materials  and  Specimen  Preparation 

Ferroelectrics  considered  for  large  displacement  actuator  applications  are  of 
three  typjes^^"^^:  i)  piezoelectrics;  ii)  ferroelectric  relaxors;  iii)  antiferroelectrics.  The 
first  category  has  been  most  popular  in  practice  and  therefore  is  emphasized  in  the 
preliminary  study.  Poled  piezoelectrics,  PZTs,  were  obtained  from  a  commercial 
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sourcet:  PZT-A  (Zr/Ti  =  52/48)  with  hard  dopants  (Ba  and  Fe)  and  PZT-B 
(Zr/Ti  =  53/47)  with  soft  dopants  (Nb).  The  relaxor  material  and  the 
antiferroelectrics  were  examined  in  a  companion  study. 

Short  compression  cuboids  with  nominal  dimensions  6  x  8  x  12  mm  were  cut 
from  bulk  ferroelectrics  and  all  non-electroded  faces  groimd,  using  a  resin-bond 
diamond  wheel.  In  some  cases,  one  of  the  side  faces  was  polished  to  facilitate 
observations  during  testing.  The  electrodes  were  connected  to  a  large  capacitor  for 
pxilarization  measurement. 


Fig.l.  Compression  loading  fixture  with  associated  measurement  circuits. 

2.2  Testing  and  Measurements 

Uniaxial  compression  testing  of  brittle  materials  is  notoriously  difficult,  given 
the  requirement  that  the  load  be  applied  at  the  geometric  center  to  eliminate 
bending.  Nearly  perfect  alignment  is  essential  in  order  to  characterize  the  non-linear 
deformation.  Consequently,  a  special  loading  fixttue  (Fig.  1)  has  been  designed  to 
facilitate  the  alignment  procedure  and  also  to  provide  a  means  for  recording  the 
electric  displacement.  Two  spherical  cones  were  incorporated  with  two  flat  alumina 
blocks  used  to  insulate  the  electrodes.  Specimens  were  carefully  emplaced  at  the 
center  between  the  alumina  blocks.  The  fixture  was  then  incorporated  into  a 
servohydraulic  machine  for  testing.  Final  adjustments  needed  to  minimize  bending 
were  made  through  in  situ  measurement  of  the  strain  differential  upon  exposing  the 

t  Channel  Industries,  Inc.,  Goleta,  California. 
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specimen  to  loading/ unloading  cycles,  within  a  very  small  stress  range  (<  5  MPa). 
The  specimen  was  connected  in  a  closed  loop  with  a  large  capacitor,  having 
capadtarce  five  orders  of  magnitude  larger  than  the  spedmen  (Fig.  1).  In  this  case, 
the  voltage  change  across  the  capadtor,  measured  by  a  voltage  follower,  provided 
an  indication  of  the  charge  transfer  from  the  specimen  to  the  capadtor. t 
Consequently,  the  polarization  change,  AP,  could  be  obtained  as 

AP  =  CqV  (1) 

where  Co  is  the  capacitance  of  the  shunt  capadto~  ind  V  is  the  voltage  measured 
across  the  capacitor.  The  convention  adopted  here  is  that  depolarization  is  indicated 
by  positive  AP. 

Other  quantities  measured  simultaneously  with  the  voltage  include  the  load 
instantaneously  carried  by  the  spedmen  obtained  from  the  load  cell  output  of  the 
testing  machine.  Also,  bi-directional  strain  gauges  were  emplaced  on  two  opposing 
side  faces  to  measure  the  deformation  along  both  the  loading  and  transverse 
directions,  as  well  as  monitor  possible  bending  caused  by  misalignment.  The  strain 

component  along  the  loading  direction  is  designated,  For  spedmens  tested  with 
the  loading  perpendicular  to  the  polar  axis,  the  strain  along  that  axis  is  designated, 

and  the  other  normal  strain,  Eyy.  When  the  loading  is  parallel  to  polar  axis, 
e»v  =  Evv  Positive  quantities  for  stress,  and  strain,  denote  compression,  unless 

XX  yy 

otherwise  indicated. 

3.  RESULTS 

The  deformation  behavior  found  for  hard  and  soft  PZTs  subject  to  stress 
applied  parallel  to  the  poling  direction  (Fig.  2)  indicate  that  both  materials  are  non¬ 
linear  when  the  stress  exceeds  -  20  MPa,  and  exhibit  similar  permanent 
deformations  (-  0.5%)  after  exposure  to  stresses  above  500  MPa.  The  major 
difference  between  these  materials  occurs  upon  unloading.  In  the  "soft”  material 
(PZT-B),  the  non-linear  deformation  is  nonrecoverable,  at  stress  levels  as  low  as 
30  MPa  (Fig.  2a).  Conversely,  the  “hard"  material,  (PZT-A),  recovers  most  of  the  non¬ 
linear  deformation,  with  considerable  hysteresis  (Fig.  2b);  a  phenomenon  usually 
referred  to  as  anelastic  deformation.  However,  the  capacity  for  anelastic  strain  is 
destroyed  at  stress  >  300  MPa. 

The  strain  measurements  indicate  that  the  volume  strains  (Ezz  £xx  +  Eyy) 
essentially  linear  for  both  materials  (Fig.  3).  Consequently,  the  non-linear 
deformation  is  related  solely  to  the  deviatoric  components  of  stress,  CT',  and  strain,  £' 
(Fig.  4).  The  instantaneous  Poison's  ratio  obtains  a  value  characteristic  of 
incompressible  plastic  flow  (0.5)  as  the  nonlinear  deformation  proceeds.  The 
deviatoric  stress/strain  curves  have  a  form  (Fig.  5)  characteristic  of  that  found  for 
phase  transformations,  consisting  of  i)  an  initial  hardening  regime  wherein  the 
hardening-rate  diminishes  as  the  stress  increases,  followed  by  ii)  a  regime  of 
accelerated  hardening,  as  a  saturation  plastic  strain  is  approached.  Consequently,  the 

stresses  O'  and  plastic  strain,  E'p,  can  be  described  by  a  bifunctional  power  law 
expression. 


f  It  is  emphasized  that  the  change  in  voltage  is  minimal. 
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Fig.  2.  Stress /strain  curves  from  hard  and  soft  PZTs  subject  to  stress  along 

polar  axis. 

o'/o;  =  (e;/e;)*  +  (e;/E;f  (2) 

where  a  and  b  are  empirical  power  law  coefficients  (a  <  1  and  b  >  1),  is  a 

reference  stress  and  €  'g  a  reference  strain.  The  magnitudes  of  the  coefficients  for  the 
hard  and  soft  materiiils  are  summarized  in  Table  I. 


TABLE  I  Constants  in  Deviatoric  Constitutive  Law  (Eqn.  2) 


a 

b 

(Hard)  PZT-A 

0.54 

9 

100 

0.62 

(Soft)  PZT-B 

0.7 

12 

75 

0.65 
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Fig.  3.  Change  in  volume  strain  with  axial  stress  for  hard  and  soft  PZT's. 


Fig.  4.  The  deviator  plastic  strain,  as  it  relates  to  deviatoric  stress,  is  solely 
responsible  to  the  non-linear  deformation. 
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Fig.  5.  Electrical  depolarization,  AP,  for  hard  and  soft  PZTs  as  a  function  of 
the  applied  stress. 

For  the  both  materials,  the  charge  on  the  capacitor  changes  during  loading, 
expressed  as  a  depolarization,  AP.  The  depolarization  curves  for  each  material 
(Fig.  5)  resemble  the  corresponding  stress/strain  curves  (Fig.  2),  suggesting  that  the 
sources  of  the  non-linear  deformation  are  also  responsible  for  the  polarization 
changes.  The  total  depolarization  measured  after  a  stress  exposure  of  -  500  MPa  is 
-  80%  of  the  remnant  polarization  induced  upon  px>ling. 

The  stress/strain  curves  obtained  when  the  stress  is  applied  normal  to  the  poling 
direction  are  qualitatively  the  same  (Fig.  6)  but  exhibit  more  limited  non-linearity. 
For  example,  the  permanent  strain  in  the  soft  material  (-  0.3%)  is  about  half  that 

obtained  upon  loading  along  the  poling  axis.  However,  the  transverse  strains  Exx 

and  Eyy  (Fig.  7)  have  opposing  non-linearities,  as  a  result  of  the  anisotropy  caused  by 
poling.  The  depolarization  curves  exhibit  features  (Fig.  8)  dramatically  different 
from  those  evident  upon  loading  along  the  poling  axis.  Initially,  the  polarization 
increases  as  the  strain  increases,  but  when  the  strain  exceeds  ~  0.1%,  an  inversion 
occurs  and  the  polarization  begins  to  decrease  with  continued  straining. 
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Fig.  6.  Stress/ strain  curves  for  hard  and  soft  PZTs  when  the  stress  is  applied 

normal  to  the  polar  axis. 

4.  NON-LINEAR  DEFORMATION 

The  volume  strain  in  all  specimens  appears  to  be  linear  and  non-hysteretic 
(unless  microcracking  occurs).  Consequently,  when  non-linearity  arises,  it  is  only 
necessary  to  understand  and  model  the  deviatoric  deformation:  as  in  metal  plasticity. 
The  incremental  deviatoric  strain  u]X)n  loading  then  relates  to  the  corresponding 
deviatoric  stress,  through  a  constitutive  law  (Eqn.  2). 

The  mechanism  responsible  for  the  non-linear  deformation  in  the  piezoelectric 
materials  involves  90®  domain  switching,  having  the  features  schematically 
illustrated  in  Fig.  9  The  as-processed  polycrystalline  material  has  a  nominally 
random  set  of  orientations  in  the  polarization  vector,  Pj.  An  applied  electric  field 
induces  both  180®  and  90®  domain  switching,  leading  to  a  preferred  alignment  of  P, 
in  the  direction  of  the  field.  Residual  stress  also  develops  in  association  with  90° 
domain  switching.  Poling  of  PZTs  is  conducted  with  a  field,  E^.  Subsequent  stressing 
may  be  along  either  z  or  x.  For  domain  B  (Fig.  9),  oriented  along  z,  initially 
Px  =  Py  =  0  and  Pz  =  Pq.  This  type  of  domain  is  most  likely  to  switch  during 
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Fig.  7.  Transverse  strain  characteristics  for  both  PZTs  as  a  function  of  applied 
normal  stress. 

poling,t  above  a  critical  field,  to  an  orientation  Py  =  =  0  and  P*  =  Pq.  The  new 

orientation  is  stable,  subject  to  subsequent  mechanical  stress  Ozz,  but  depoling  is 
possible  when  a  <Jxx  stress  is  applied.  If  poling  did  not  cause  the  domain  to  switch 

(because  the  field  is  too  small),  switching  may  be  activated  by  a  subsequent 
stress,  whereupon  the  domain  would  realign  along  either  P*  or  Py.  For  domain  C,  the 
situation  is  similar  and  an  electric  field  Ex  is  again  likely  to  induce  Pv  — »  Px 
switching.  Domain  D  is  most  illustrative  of  the  interplay  between  the  electric  field 
and  stress.  No  domain  switching  can  occur  subject  to  Ex,  because  there  is  no  net 
change  in  Px  upon  switching  from  P,  to  P2:  moreover,  P*,  decreases  if  switching 

occurs  from  Pj  to  P3.  However,  upon  Czz  stressing,  switching  from  Pj  to  Pj  is 
expected  when  the  normal  stress  (a^  =  Ozz/2)  reaches  a  critical  value.  This  event 


f  In  some  cases,  90°  domain  switching  upon  poling  may  be  inhibited  by  the  strain  energy  increase 
associated  with  the  residual  strain. 
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produces  a  positive  strain  i*'  6  yy  I’ut  negative  strain  in  and  A  related  mode  of 
switching  may  occur  upon  <Jxx  stressing. 

The  basic  insight  obtained  from  Fig.  9  allows  the  major  phenomena  found  by 
experiment  to  be  rationalized.  For  the  pded  PZT  materials,  tet^  with  the  mechanical 
loading  axis  coincident  with  the  poling  direction,  the  mechanical  field  causes 
switching  in  a  sense  opposite  to  that  caused  by  poling  resulting  in  systematic 
depoling.  Once  the  depoling  is  complete,  further  switching  does  not  occur. 
However,  this  mechanically  depoled  state  must  be  contrasted  and  distinguished 
from  the  origiiuLl  (unpoled)  state,  in  that  the  random  polarization  vectors  in  the 
virgin  material  are  now  in  a  plane  normal  to  the  depoling  stress.  "Complete 
depoling"  results  in  the  saturation  behavior  found  at  plastic  strains  >  0.5%  (Fig.  2). 
Thereafter,  the  material  behaves  as  an  elastic  body  until  microcracking  commences 
prior  to  mechanical  failure. 


Fig.  8.  Polarization  changes  in  the  polar  axis  caused  by  stress  applied  normal 
to  the  polar  axis. 
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Fig.  9.  Schematic  of  deformation  by  domain  switching. 

Loading  normal  to  the  poling  axis  (usually  thought  to  be  linear)  causes  a  non¬ 
linear  response  best  illustrated  by  orientation  D  in  Fig.  9.  Some  of  the  90°  domains 
are  unfavorably  oriented  for  switching  upon  poling,  but  are  in  a  preferred 
orientation  for  switching  upon  loading.  The  important  consequence  is  the 

relationship  found  between  the  polarization  change,  AP,  and  the  strain  (Fig.  8). 
The  initial  portion  of  these  curves  is  the  normal,  linear  piezoelectric  effect.  However, 
type  D  domain  switching  at  larger  stresses  causes  depolarization  (-  APx)  and  an 
expansion  along  the  x-axis,  leading  to  the  observed  inversion  in  the  polarization 
curve. 

The  parameters  that  characterize  the  non-linear  deviatoric  stress/strain  (Eqn.  2), 
must  depend  on  factors  such  as  the  domain  wall  mobility,  the  residual  stresses,  etc. 
More  detailed  analysis  is  needed  to  develop  a  general  yield  criterion  and  a  flow  law. 
For  this  purpose,  experiments  involving  multiaxial  stresses  are  necessary. 
Furthermore,  the  differences  found  between  the  "hard"  and  "soft"  PZTs  should  be 
consistent  with  the  yield  criterion  and  flow  law. 

5.  CONCLUSION 

Mechanical  non-linearity  found  in  piezoelectric  ceramics  has  been  attributed  to 
ferroelastic  domain  switching,  with  a  yield  criterion  controlled  by  a  deviatoric 
effective  stress.  The  flow  rule  obeys  a  power  law  formula  at  small  strains,  similar  to 
metal  plasticity.  At  larger  strain,  accelerated  hardening  occurs  and  the  inelastic 
strain  approaches  saturation.  Consequently,  a  bifunctional  power  law  is  proposed  to 
model  the  overall  nonlinear  response. 


( 


Strong  interconnections  have  been  established  between  ferroelectric  and 
ferroelastic  deformation.  Electrically  "hard"  PZT  can  sustain  appreciable  mechanical 
loads  without  signiftcant  permanent  deformation,  whereas  "soft"  PZT  is  susceptible  to 
permanent  deformation  at  small  stresses.  This  difference  arises  even  though  the 
"yield"  stresses  are  similar  for  the  two  materials.  Consequently,  for  actuator 
applications,  hard  PZT  can  be  used  at  higher  stress  levels,  without  dep>oling. 
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ABSTRACT 


Singular  stress  and  electric  fields  around  a  crack  are  studied  for  a  linear  elastic, 
electrostrictive  body  by  using  an  energy  conservation  integral  and  by  solving  the  eigein 
functions  to  the  equilibrium  equations.  When  the  dielectric  permeativity  of  the 
materials  is  very  large  (>106  - 10^  eq  ),  the  stress  and  electric  field  decouple  and  exhibit 
the  classical  square-root  singularity  (r  *V2).  However,  when  the  electrostrictive 
deformation  becomes  comparible  with  the  elastic  strain  at  the  crack  tip,  the  asymptotic 
solution  obtained  form  this  study  suggests  that  the  stress  and  electric  field  exhibit 
different  types  of  sigularity:  with  stress  exhibiting  square-root  singularity  ( r  )  and 
the  electric  field  exhibiting  a-quarter-root  singularity  ( r  ).  Complete  angular  field 
distributions  are  obtained  numerically  for  the  asymptotic  solutions.  Energy  release 
rates  are  also  derived  to  correlate  the  intensity  factors. 

Prelimenary  comparison  of  the  solution  with  recent  experiments  on  crack  extension 
in  an  electrostrictive  material  is  presented.  These  materials  are  intended  for  advanced 
actuator  applications.  Implications  for  mechanical  reliability  of  multilayered 
ferroelectric  actuators  are  discussed. 
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NOMENCLATURE 


Di 

Ei 

Ei(e) 

eij  (0) 

G 

/(v,m) 

J 

K 

Ke 

Kie 

m 

Pi 

Qii  'Qi2 

(r,0) 

S 

Si,S2 

Ui 

(x,y) 

Y 

yi^-ye 


Crack  length 

Electric  displacement  vector  ( electric  induction ) 

Electric  field  vector 

Angular  function  of  strain  around  a  crack  tip 
Strain  tensor  ( small  strain  theory) 

Angular  function  of  strain  around  a  crack  tip 
Shear  modulus 
Numerical  fimction 
J-integral  ( energy  release  rate) 

Jacobian  for  planar  approximation  of  surfaces 
Dielectric  constant 
Electric  intensity  factor 

Stress  intensity  factor  induced  by  electromechanical  coupling 
Electrostriction  ratio;  m=S2/S| 

Electric  polarization  vector 

Electrostrictive  constcmts 

Polar  coordinate  with  crack  tip  as  origin 

Residuals  at  the  left  boundary  condition  for  shooting  method 

Equivalent  electrostrictive  coefficient 

Electrostrictive  coefficient 

Mechanical  displacement 

Cartisian  coordinates 

Young's  modulus 

variables  for  numerical  solutions 

Kronecker-delta 
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e  Electric  permitivity  (e=  K  Eq  ) 

Eq  Electric  permitivity  of  vacume  (£(,  =  0.88  X  10'^^ ) 

'E  Potential  energy  of  a  body 

0  Electric  potential 

Ojj  Mechanical  stress  tensor 

V  Poisson's  ratio 
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1.  INTRODUCTION 


Electrostrictive  materials  has  been  widely  employed  as  capidtor  material  where 
application  demands  large  capadtance  but  limit  space  is  available.  More  recently,  they 
are  been  designed  for  large  strain  actuators  to  be  incorporated  into  smart  materials  and 
systems^'2.  In  these  applications,  they  are  subjeded  to  large  mechanical  loading  as  well 
as  the  high  electric  driving  field.  Consequentely,  the  mechanical  and  electrical 
degradation  has  become  an  important  concern  3.  McMeeking  has  attempted  to  analyze 
the  stress  distribution  around  a  crack-like  flaw  in  a  linear  elastic,  electrostrictive  body 
using  an  energy  conservation  integral  approach^'^  solutions  are  limited  to  less 
sharp  flaws,  such  as  processing  porosities,  where  the  stress  concentration  is  less  severe 
than  those  found  around  a  sharp  crack.  The  stress  and  electrical  field  distribution 
around  a  brittle  crack  have  not  been  rinvestigated.  'ent  of  this  artical  is  to  derive 
the  solutions  in  a  rigorious  and  self-con'^isteiti  fashion,  such  that  a  plausible  fracture 
criterion  could  be  derived  and  furfhermore  be  used  to  predict  the  mechanical  failure  in 
these  materials. 


2.  CONSTITUTIVE  EQUATIONS 


We  consider  a  linear  elastic  electrostrictive  dielectric  body.  The  constitutive 
properties  can  be  described  in  two  dimensional  problem,  by  the  following  equation;^ 


-  ~  (Qii “  Qi2)PiPj  QizPkPk^ij 


J  ■  ■  (1) 

where  o  is  the  mechanical  stress,  e  is  the  inifinitesimal  strain,  P  is  the  local  electric 
palorization,  G  is  the  elastic  shear  modulus,  v  is  the  Poisson's  ratio,  and  Q12  are  the 
electrostrictive  constants  and  5ij  is  the  Kronecker's  delta  ( 5ij=l  if  i=j,  and  5ij=0  if  i?tj ). 


1  In  all  equations,  summation  over  repeated  index  is  implied  .  Index  1  and  2  refer  to  coordinates  x  and 
y  in  Cartesion  coordinate  and  r  and  0  in  polar  coordinate. 
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Here  the  term  mechanical  stress  referes  to  tractions  transmitted  by  contact  in  a  body. 
Consequently  the  Comlombic  forces  are  treated  as  body  forces  ( which  equals  to  the  flux 
of  Maxwell  stresses  in  a  rigid  dielectric  body).  If  the  material  behaves  as  a  linear 
dielectric  under  stress-free  state,  that  is; 

Pi=eEi  (2) 

where  e  is  the  free-body  dilectric  permeativity  and  £=K£fy  with  being  the  dielectric 
permeativity  of  the  vacume  cind  K  the  dielectric  constant.  Equation  (1)  can  be 
rearranged  by  making  use  of  eqns.  (2); 


Oij  =  2G|^eij  +  J  +  S^EjEj  +  S2EkEk5ij 

where  S|  and  S2  are  the  electrostrictive  coefficients,  which  can  be  expressed  as; 

51  =  -2G  (Qn  -  Q12) 

52  =  -2G(Q„  +  Q,2)^e2 

Under  the  action  of  both  the  mechanical  strain  and  electric  field,  the  electric 
displacement  is  related  to  E  by; 

Dj  =  EjjEj  =  (cSjj  —  2S^ejj  —  2S2eyj5jjj  Ej 


3.  ENERGY  DENSITY  AND  CONSERVATION  INTEGRAL 

The  mechanical  and  electric  equilibrium  conditions  provide  the  field  governing 
equations; 

Oii,i  =  0 

Du  =  0 

where  no  body  force  or  free  charges  are  present  in  the  body.  A  potential  energy 
expression  is  sought  whereby  the  conjugate  variables  (E,  D)  and  (a,  e)  are  related  by; 
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Integration  of  eqns.  (3)  and  (5)  leads  to; 


'F  =  G 


(8) 


V  2  ^  /  V  g 

■^(EiEj  +  mEj,Ej,5ij)  eij-  —  E^Ek 
McMeeking  has  derived  an  energy  conservation  integral  for  two  dimensional  problems; 
J  =  Jj,  ['Eni  -  njtjkUk,!  +  njDjEj  ds 


(9) 


In  order  that  J  is  independent  of  any  contour,  F,  enclosing  the  crack  tip,  it  is  necessary 
that  the  last  two  terms  vanish  on  the  crack  surfaces. 


3ii  dUv 


(10) 


It  is  clear  that  both  a  traction-free  conducting  crack  and  a  traction-free  impermeable 
crack  satisfy  eqn.  10.  In  the  former.  Ex  vanishes  along  crack  surface  (O  =  constant).  And 
in  the  latter,  Dy  =  0  along  the  crack  surface. 


4.  CRACK  TIP  SINGULARITIES 


Evaluating  J  integral  from  eqn.  9  over  a  circular  path  of  radius  r  from  the  crack  tip 
leads  to, 

-=  J  j^[e(r,0)]cos(0)-T(r,6)^(r,e)  +  [DrCOs(0)  +  Desiw(0)]Ex|ds 
^  (11) 

Evidently,  for  J  to  have  a  finite  value,  the  terms  in  {...)  must  exhibit  a  sigularity  as  r 

approaches  the  crack  tip.  A  carefull  scrutiny  of  these  terms  suggests  that  they  consist  of 

three  different  combination  terms: 
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^  “  Gejjeij  +  SiEiEjCij  +  EEkEk 


sG 


eijeij  + 


Ej  L  ,  Ek 


(12) 


Three  classes  of  sigularities  might  be  exhibited  by  the  strain  and  the  field  components. 
They  will  be  discussed  consecutively  in  the  following. 

Case  I: 

For  many  advanced  relaxor  ferroelectric  ceramics  (  such  as  PMN«PT ),  the  dielectric 
constant  is  on  the  order  of  10^.  Consequently,  it  is  nature  to  assume  that  the 
contribution  to  dielectric  permeability  from  the  mechanical  distortion  through  the  last 
two  terms  in  eqn.  5  can  be  ignored  compared  with  e.  However,  the  strain  produced  by 
the  electric  distortion  is  significant  when  compared  with  the  mechanical  strain.  The 
electric  force  can  be  simulated  by  a  body  force.  In  this  case,  the  e  and  E  both  possess 
conventional  sqare  root  sigularity; 

Ei~KEr-''2E,(e)  (j3, 

where  Kjg  and  Kg  determines  the  amplitude  of  the  field  singularity,  Cjj  and  Ejare 
angular  distribution  functions.  In  particular,  are  given  by  the  conventional  Mode  I 
elastic  solutions  and  Ej  coincide  with  the  angvilar  functions  in  the  Mode  IH  elastic 
solutions.  McMeeking  has  analyzed  this  situation  where  he  provided  the  following 
solutions  for  a  crack  of  length  2a  in  an  infinte  body  subject  to  an  uniform  remote  electric 


field  E: 


KjE  =  -4S(Ey)^VTW 
Ke=e;;V^ 


where  S  =  - (S.  +  2S2) 

8(1 -v)^  ^ 


(14) 
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It  must  be  noted  that  for  many  electrostrictive  ferroelectrics,  the  cx)effiecient,  S,  is 
negative  by  defination  (Eqn.  4),  and  therefore  the  stress  intensity  factor  is  positive.  In 
contrast,  for  some  common  dielectrics  such  alunnina  and  quartz,  S  is  positive  and  the 
electric  field  is  said  to  suppress  the  mechanical  cracking. 

Case  11; 


If  the  strain,  e^,  and  the  field,  Ej,  are  both  singular,  the  third  term  in  the  energy 
density  function  ( eqn.  8),  'F,  would  become  negligibly  small  as  r  goes  to  zero. 
Consequently,  we  can  make  the  first  and  the  second  term  to  exhibit  the  same  oredr  of 
singularity. 


and 

E,=o(rl) 

Substitution  of  eqns.  15  into  eqn.  8  leads  to; 

'F  =  0(r2p)  +  0(r2^-"P) 

Therfore,  we  must  have:  p=-l/2  and  q=-l/4.  Consequently,  the  following  formal 
solutions  are  sought; 


where  Uj.,  uq,  4>  are  angulcir  functions  to  be  determined  by  the  field  governing 
equations  (eqn.  6). 


(15) 


(16) 
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5.  STRESS  AND  FIELD  AT  AN  IMPERMEABLE  CRACK 


In  this  section,  an  impermeable  crack  with  type  n  singularity  is  considered.  The 
solution  can  be  expressed  by  eqn.  16.  The  strain  and  electric  field  can  be  derive  from 
these  expressions; 


8u- 

“  r  r  ae 


1 

e,e-- 


ar 


r:  ^ 

^-‘aT 

c 

®  r  ae 


1  au, 

+  — 


ae 


(21) 


The  stresses  and  electric  displacements  can  be  obtained  by  subsituting  eqn.  21  into  eqn. 
3  and  5.  Specifically,  the  following  equations  exist; 


a„  =  2G  e,^  + 


<j0e  =  2G  e00  + 


"  l-2v 

V 


^kk 


+  SiErEr  +  S2EkEk 


®kk  ■'■SiEeEe  +  S2EjjEi^ 


l-2v 

==  2Ge,fl  +  S-iErEe 

Substitution  of  eqns.  16  and  21  into  eqn.  22  leads  to; 


Orr  = 


<^00- 


"  -f  ■  ■■  U0  +  —(1  +  m)C>^  +  mO'^ 

1- 2v  "  l-2v  ®  16^  ^ 

2- v  -  2(l-v)_,  9m  =2  y*/2 

- u,  +  -^ - ^Ue  + - 0  +(l  +  m)<I) 

l-2v  '  l-2v  ®  16  ^  . 

ii:--U0  +  -<i4»' 

.  2  ®  4 

Dr  =  [4  3m)urO  -|c>'ue  +  ^Ou^  +  u;0'j 


f— 1 

1 _ 

V  r  J 

L 

rjG> 

1  r  > 

(?: 

1 - 1 

j-JG' 
1  r  / 

(22) 
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(23) 


D,  =  (2+ 3mK®'+ |u;a.  +  2(1  +  m)a^ *'  -  |ue« 

where  m=S2/Si.  From  the  equalibrium  condition  (eqn.  6),  we  have  explidtely; 

■  <^rr  ~  <^ee  .  ^  _  q 

dr  r  r  ae 

+ ZEifi + =  0 

dr  r  r  d0 

I  I  ^  ^De  _  Q 

dr  r  r  d0  (24) 

Finally,  the  governing  equations  for  the  angular  functions  are  obtained,  by  substituting 
eqn  23  into  eqn  24,  as; 


3  1-v  .  5-8v  .,  .  9('l-m)^2  l  +  2m^,2  „ 

Ur - U. - ; - tUq  + - <D - O  +  — <1)0  =  0 

"  21-2V  "  2(l-2v)  ®  32  4  4 


31-2V.  .  7-8v  9a  +  m)(l-2v)^^,  a  +  m)(l-2v)^,^„  ^ 
®  8  1-V  ®  4(l-v)  "  16(l-v)  1-v 

160"[(2  +  3m)Ur  +  2(1  +  m)u^]  +  80'[4(1  +  m)u^'-  Uq  +  6(1  +  m)u;] 

+  <i>[l2ur  +  3(1  +  3m)ur  -  6(1  -  m)ue]  =  0 

To  fadlitate  the  numericed  scheme,  eqn.  25  can  be  rearranged  into  a  matrix  form; 

10  -O  1 

4  ul' 

(l  +  mXl-2v)-, 

0  1  - '-<!>'  Ue  = 

1-v  . 

120  32(1  +  m)0'  16[(2  +  3m)ur  +  2(1  +  m)u^] 


3  1-v  _  5-8v  _,  9(1 -m) -2  l  +  2m -  ,2 

21-2V  "  2(l-2v)  *  32  4 

31-2V.  7-8v  .,  9(l  +  m)(l-2v)^, 

8  1-v  4(l-v)  ^  16(l-v) 

80'[ug  -  6(1  +  m)u']  -  0[3(1  +  3m)u^  -  6(1  -  m)ue] 

Now  consider  the  boundary  conditions.  At  0=0,  synunetry  condition  requires; 


(25a) 
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and 


a. 

b. 

c. 

d. 

e. 


0^  =  0 


30i 


de 

3a, 


fla  =  o 


36 


rr  « 


=  0 


Dr  =  0 


3De_ 


=  0 


36  (26) 

These  conditions  can  be  satisfied,  through  symmetry  argument,  by  requiring;  at  6  =  0 

U0  =  O 

u;  =  o 


0  =  0 


(27) 


In  addition,  for  a  traction-free,  impermeable  crack,  we  must  have  on  the  crack  surface 
(at  6  =  7C ); 


a.  000  =  0 

b.  0^  =  0 

c.  D0  =  O 


(28) 


The  assiuned  forms  of  the  solution  (eqn.  16)  warrante  path-independence  of  the  J- 

integral,  validating  eqn.  10  for  all  solutions.  Consequently,  only  two  of  the  eqns.  28  need 

to  be  enfoced.  We  choose  (a)  and  (b),  i.e.  at  6  =  k  : 

u;--U0  +  -OO'  =  O 
"  2^  4 


2-v  _  2(l-v).,  9m  =2  /i  \^f2  o 

- u,  +  -^ - + - O  +(l  +  m)0  =0 

l-2v  '  l-2v  ^  16  ^  ' 


(29) 


Additionally,  because  the  equations  are  homogeneous,  a  normalization  is  introduced 
such  that  at  6  =  0; 

^  =  (30) 

Eqnuations  (27),  (29)  and  (30)  provide  the  complete  set  of  boundary  conditions. 
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6.  NUMERICAL  SOLUTIONS 


Numerical  solution  to  eqn.  (25)  with  bondary  conditions  (eqns.  28, 31  and  32)  is 
pjerformed  by  means  of  a  secant  method  for  solving  nonlinear  systems  of  equations. 
Regard  the  problem  as  an  initial  value  formulation  subject  to  initial  condition  at  0  =  0, 

Ur  = 

(31 

where  w^'^  and  W2^  are  two  initial  guesses.  Deviation  from  the  BC  eqn.  30  is  measured 

by  the  two  remainders  defined  through; 

R^/^  =  u;--Ue  +  -0«i>' 

2  4  l0_,r 

rH>  =  Izl.  u  +  uj  +  —  +  (1  +  m)4>'^ 

'  l-2v  1-2V  *  16  '  '  I,..  ,32 

For  numerical  implementations,  the  following  symbles  have  also  been  used; 

■y,l  pr' 

Yi  Ue 

4) 


Ys  Ue 


By  means  of  a  shooting  scheme,  the  corresponding  initial  conditions  are; 


yi 

Y2 

Y3 

Ya 

Ys 

-y6.e=o 


0.73275 


0.92785 


It  is  noted  that  the  convergent  solution  renders  the  system  of  equations  (25)  singular. 
Consequently,  a  Taylor  series  expamsion  form  was  impletemented  to  obtain  the  correct 
solution.  The  graphical  representation  of  the  results  are  displayed  in  Fig.  2. 
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With  the  angular  distribution  functions  displayed  in  Fig.  2,  the  J-integral  cab  be 
explidtely  evaluated,  leading  to; 

J  =  K2+...  (35) 

It  is  further  noted  that  the  conventional  crack  tip  elastic  stress  solution  is  obtained  by 
using  the  following  conditions; 


>1' 

'2' 

0 

ys 

0 

74 

0 

75 

-1 

.76. 

0=0 

0 

(36) 


8  IMPLICATION  AND  CONCLUDING  REMARKS 

Crack  tip  field  solutions  have  been  obtained  for  an  impermeable  crack  in  a  linear 
elastic  electrostrictive  ferroelectric  ceramic  body  subject  to  electric  boundary  conditions 
at  remote  locations.  In  the  asymptotic  form  as  the  crack  tip  is  approached,  the  stresses 
exhibit  conventioncd  sqare-root  singularity  ( «=  r^/2)  while  the  electric  fields  exhibit  one- 
quarter-root  singularity  ( «  Energy  release  rate  in  the  conventional  sense  has  also 
been  derived  and  related  to  the  stress  and  field  intensity  factors.  It  is  suggested  that  the 
failure  can  occur  due  to  filed  induced  stress  intensifications. 

Preliminary  experiments  have  been  conducted  using  a  relaxor  ferroelectric  ceramic 
body.  The  observations  indicated  that  the  crack  had  indeed  propagated  under  the 
imposed  electric  field  alone.  Using  a  narrow  beam  specimen,  it  is  estimated  that  the 

field  concentration  factor  at  failure  approaches _ .  In  comparison,  when  a  sp)ecimen 

was  subjected  to  mechanical  loading,  the  fracture  toughness  measured  at  Kjc  =  0.9 
MPa  Vm. 
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Table  I 


Material  Properties  for  .9PMN».1PT 


Young's  modulus,  Y  (GPa) 

125 

Shear  Modulus,  G  (GPa) 

46 

Poisson's  ratio,  v 

0.35 

Dielectric  constant,  K 

5.25  X  104 

Electrostrictive  constants. 

Q„  (m4/C2) 

0.0137 

Qi2  (mVC2) 

-0.00514 

Electrostrictive  coefficients. 

Si  (mVC2) 

-3.7x10-3 

-1.7x10-3 

S(,r  /  C2) 

-4.5310-5 

Fracture  toughness,  Jic  ( 1/1^12 ) 

8 
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Figure  Captions 

Fig.  1  Crack  geometry  in  a  dielectric  material 

Fig.  2  Angular  distribution  functions  for  a  crack  in  a  linear  elastic  electrostrictive 
material. 

Fig.  3  Crack  extension  observed  in  PMN  •PT  under  alternate  electric  field(~3MV  / m). 
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